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ABSTRACT 


In  the  past,  contagious  distributions  have  been  successfully 
applied  in  Bacteriology,  Entomology  and  Accident  Statistics. 
This  thesis  applies  the  notion  of  contagious  distributions 
in  the  inventory  control  of  new  products  and  seasonal  or 
style  goods,  which  have  an  underlying  "true  contagion"  for 
their  demands,  viz,  the  influence  of  past  demands  on  future 
occurrence  of  demands. 

A  contagious  distribution  is  derived  by  assuming  a  modified 
Poificon  process  where  the  demand  rate  at  any  instant  of  time 
depends  on  the  past  demands  prior  to  that  instant.  A 
discussion  on  estimation  of  the  various  parameters  of  the 
contagious  distribution  is  also  included.  Using  this 
contagious  distribution,  a  multi-period  inventory  model  is 
discussed  for  new  product  lines  with  a  "fixed  periodic 
review  policy."  An  optimal  s  -  S  order  policy  is  derived 
as  a  function  of  the  initial  stock  level  and  the  review 
period.  Seasonal  or  style  goods  are  treated  as  single¬ 
period  inventory  problems  with  contagious  demands.  An 
algorithm  is  developed  to  compute  the  optimal  order  policy 
and  the  optimal  length  of  the  period. 
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INTRODUCTION 


In  literature,  successful  applications  of  contagious  distributions  are 
found  in  Bacteriology,  Entomology  and  Accident  Statistics.  In  this  thesis, 
an  attempt  has  been  made  to  apply  the  notion  of  contagious  distributions  in 
the  inventory  control  of  consumer  products  which  have  an  underlying  "true 
contagion"  for  their  demands,  particularly  when  the  product  is  new  or  is  a 
seasonal  or  style  good.  "True  contagion"  means  the  probability  of  a  favor¬ 
able  event  depends  on  the  occurrence  of  previous  favorable  events.  Thus, 
unlike  the  classical  Poisson  demand,  the  contagious  demand  in  two  nonover¬ 
lapping  time  intervals  are  dependent  as  occurrence  of  a  demand  influences 
further  occurrence  of  future  demands. 

Chapter  1  discusses  a  "contagious"  demand  model  for  an  inventory  system. 
Starting  with  literature  review  on  contagious  distributions,  the  various 
applications  of  our  contagious  demand  model  are  discussed.  It  is  to  be  noted 
that  the  contagious  behavior  is  only  transient  and  once  the  product 
stabilizes  in  the  market,  the  contagion  effect  vanishes.  By  assuming  a 
modified  Poisson  process  where  the  demand  rate  at  any  instant  of  time  depends 
on  the  number  of  past  demands,  the  contagious  demand  distribution  is  derived. 
Since  the  life  of  new  product  lines  are  longer,  compared  to  seasonal  goods 
like  clothings,  both  commodities  are  treated  separately.  Seasonal  goods  arc 
considered  as  single  period  inventory  problems  as  the  period  is  small  while 
new  product  lines  are  considered  as  multi-period  inventory  problems  composed 
of  a  number  of  single  periods.  The  contagious  demand  distributions  for 
succeeding  periods  for  new  product  lines  are  also  derived. 

Chapter  2  discusses  the  case  of  new  product  lines  exclusively.  Here,  a 
fixed  periodic  review  policy  is  followed  after  fixing  the  review  periods 
institutionally  and  solving  the  multi-period  problem  by  successively  solving 
a  single  period  problem,  knowing  the  review  period.  An  s  -  S  policy  is 
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derived  as  an  optimal  order  policy  in  this  case.  An  algorithm  to  compute  the 
optimal  order  level  for  a  given  review  period  is  also  developed.  The  optimal 
ordering  policy  as  a  function  of  the  length  of  the  review  period  is  also 
plotted. 

Chapter  3  exclusively  discusses  the  seasonal  or  style  goods  which  are 
considered  as  single  period  inventory  problems.  Besides  a  fixed  period- 
length  policy,  an  optimal  period-length  policy  is  also  discussed.  For  the 
latter  case,  an  algorithm  to  compute  the  optimal  review  period  is  developed 
after  computing  the  optimal  order  level  for  a  given  review  period.  To  the 
author's  knowledge,  this  is  the  first  time  an  optimization  over  review  period 
is  successfully  carried  out  for  an  inventory  policy  under  stochastic  demand. 

Chapter  4  discusses  two  methods  of  statistical  estimation  of  the  param¬ 
eters  of  the  contagious  demand  distribution.  Here  it  is  also  shown  that  the 
maximum  likelihood  estimator  is  asymptotically  unbiassed  and  efficient. 

Appendices  on  solution  of  ordinary  differential  equations,  properties  of 
Beta  and  Gamma  functions  which  appear  in  our  analysis  are  also  included  as 


reference. 


CHAPTER  1 


A  "CONTAGIOUS"  DEMAND  MODEL 

1.1  Literature  Review  on  Contagious  Distributions 

In  a  recently  published  volume  [10],  of  considerable  interest,  covering 
the  Proceedings  of  the  International  Symposium  on  Classical  and  Contagious 
Discrete  Distributions,  we  observe  the  successful  applications  of  the 
notion  of  contagious  distributions,  particularly  to  biological  populations, 
accident  statistics,  contagious  diseases,  and  psychological  data.  Actually, 
the  interest  on  contagious  distributions  dates  back  to  1920  when  Greenwood 
and  Yule  [6]  developed  a  very  general  scheme  for  contagious  events  where 
the  occurrence  of  each  event  increases  (or  decreases)  the  probability  of 
further  events.  But  due  to  the  very  generality  of  their  model,  their 
formulas  become  too  complex  for  practical  applications. 

J.  Neyman  [9]  developed  and  applied  successfully  three  types  of 
contagious  distributions  in  Entomology  for  distribution  of  larvae  in 
experimental  plots  and  bacteriology.  In  a  follow  up  paper  by  W.  Feller  [5], 
it  was  pointed  out  that  there  are  two  kinds  of  contagion  as  "true  contagion" 
and  "apparent  contagion"  and  Neyman' s  [9]  contagious  distributions  are  of 
the  latter  type. 

Apparent  contagion  is  the  result  of  inhomogeneity  arising  from 
distributions  on  the  parameters  involved  in  a  population.  Thus,  the 
compound  Poisson  distribution  is  an  example  of  apparent  contagion  as 
developed  by  Greenwood  and  Yule  [6]  and  applied  successfully  in  accident 
proneness. 

In  the  case  of  true  contagion,  the  probability  of  a  "favorable"  event 
depends  on  the  occurrence  of  previous  favorable  events.  Ironically  enough, 
assuming  true  contagion,  Eggenberger  and  Polya  [3]  arrived  at  the  same 
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distribution  as  obtained  by  Greenwood  and  Yule  [6].  The  Greenwood-Yule- 
PoJ ya-Eggenborgtr  distribution,  which  is  a  negative  binomial  distribution, 
has  found  many  applications  in  contagious  diseases,  sickness  and  accident 
statistics.  Curland  [7]  discusses  a  survey  of  the  applications  of  the 
negative  binomial  and  other  contagious  distributions  with  special  reference 
to  some  medical  data. 

1.2  Occurrence  of  Contagion-Demand  in  Practice 

So  far  in  the  literature  on  contagious  distributions,  no  attempt  has 
been  made  to  apply  the  notion  of  contagious  distributions  in  the  inventory 
control  of  consumer  products  which  have  an  underlying  true  contagion  for 
their  demands,  particularly  when  the  product  is  new  or  is  a  "style"  good 
(which  changes  its  style  periodically).  The  classical  demand  distribution 
used  in  inventory  control,  like  the  Poisson  distribution,  will  not  reflect 
the  true  behavior  of  the  contagious  demand.  It  is  well  known  that  the  simple 
Poisson  distribution  describes  mutually  independent  events;  in  other  words, 
with  a  Poisson  distribution  the  number  of  events  in  two  nonoverlapping  time 
intervals  are  uncorrelated  and  the  occurrence  of  an  event  has  no  influence 
on  the  probability  of  occurrence  of  further  events. 

If  we  study  the  demand  for  the  new  products  and  "style"  goods,  we  will 
note  that  besides  the  constant  demand  for  the  products  (which  is  mainly  due 
to  advertisement),  a  contagious  demand  also  occurs  due  to  the  customers  who 
have  used  the  product  and  recommended  it  to  their  friends  or  to  other 
sources  of  consumer  awareness.  So  until  the  new  product  is  stabilized  in 
the  market,  there  is  a  "contagious"  demand  during  the  initial  or  "transient" 
stage,  besides  the  constant  demand.  Unlike  Poisson  demand,  the  contagious 
demand  in  two  nonoverlapping  time  intervals  are  dependent  as  occurrence  of  a 
demand  influences  further  occurrence  of  future  demands.  Given  that  a  demand 
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has  occurred  for  a  new  product  from  a  particular  area,  we  expect  to  find 
more  demands  to  come  from  the  same  area  due  to  the  influence  of  the  past 
occurred  demand.  In  other  words,  our  demand  model  has  the  underlying 
assumption  that  past  demands  have  an  influence  on  the  occurrence  of  future 
demands. 

1.3  Applications  of  Contagion-Demand 

A  number  of  examples  of  products  can  be  thought  of,  which  will  follow 
a  contagious  law  for  their  demand.  For  example,  acquiring  a  Princess 
telephone  in  your  home  will  tempt  your  neighbor  to  do  the  same.  Other 
examples  may  be,  demands  for  new  cereals,  new  freeze-dried  coffee,  new  books, 
cars  (style  goods),  Christmas  trees  and  practically  all  consumer  products  in 
daily  use.  The  notion  of  contagion  can  even  be  extended  to  nonconsumer 
products  which  exhibit  a  contagion  pattern  like  research  reports.  In  other 
words,  a  research  organization  can  do  a  better  control  of  their  inventory  of 
reports  assuming  a  contagious  demand  distribution  for  their  reports.  Also, 
the  sequence  of  published  research  papers  in  a  particular  subject,  will  tend 
to  follow  a  contagious  law.  This  should  be  so  as  readers  after  seeing  a 
research  paper  tend  to  work  more  in  the  same  area  which  results  in  more 
papers.  This  may  be  useful  to  a  librarian  or  publisher  of  a  technical 
journal.  Another  classical  application  is  the  efficient  usage  of  hospital 
beds  for  patients  with  contagious  disease. 

1.4  A  Contagious  Demand  Model 

The  demand  model  assumes  that  every  occurrence  of  demand  produces  one 
unit  of  demand  implying  the  demand  has  a  discrete  probability  distribution. 

It  was  discussed  in  Sections (1. 2)  and  (1.3)  how  the  contagious  demand  arises 
in  practice  and  its  various  applications.  To  incorporate  this  idea  in  our 


demand  model,  wc  follow  Cox  and  Miller  [1)  by  assuming  a  constant  demand 
rate  and  superimposed  on  that  a  unit  contagious  demand  rate  to  reflect  the 
influence  of  each  past  demand.  So  at  any  point  in  time  "t,"  the  total 
contagious  demand  rate  will  be  equal  to  the  unit  contagious  demand  rate 
times  the  past  demands  before  "t."  So  the  rate  at  which  demand  will  occur 
at  any  time  "t,"  will  be  the  sum  of  constant  demand  rate  and  the  total 
contagious  demand  rate.  Thus,  we  also  introduce,  explicitly,  time  as  a 
parameter  In  our  demand  model.  Though,  the  final  contagious  demand 
distribution,  like  many  other  contagious  distributions,  is  negative  binomial 
in  some  sense,  a  derivation  of  the  distribution  will  be  carried  out  in  the 
analysis  to  follow,  to  introduce  "time"  explicitly  since  "time"  never  appears 
as  an  explicit  parameter  in  classical  negative  binomial  distributions.  (The 
reader  is  referred  to  Feller  [4]  for  a  discussion  on  classical  negative 
binomial  distribution.) 

Some  caution  must  be  applied  as  to  the  duration  and  magnitude  of 
contagious  demand.  One  cannot  expect  the  same  amount  of  contagious  demand 
throughout  the  life  of  a  product.  As  a  matter  of  fact,  the  contagious  demand 
will  be  very  high  when  the  product  has  just  been  introduced  in  the  market  and 
will  follow  a  decay  law  such  that  its  effect  vanishes  after  a  time  the 
product  is  stabilized  in  the  market.  So  the  unit  contagion  rate,  which 
reflects  the  increase  in  the  constant  demand  rate,  is  a  decreasing  function 
of  time. 

1.5  Periodic  Review  Policy 

There  are  two  types  of  phenomenal  situations  that  arise  in  practice. 


One  is  the  consideration  of  new  product  lines  like  new  cereals,  while  the 
other  is  the  consideration  of  seasonal  or  style  goods  like  automobiles  and 
clothes.  If  we  assume  that  the  contagion  rate  varies  at  each  instant  of 
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time,  this  will  complicate  our  analysis  and  results  considerably.  So,  a 
simplifying  approximation  is  made  about  the  contagion  rate  for  the  two  class 
of  goods  which  exhibit  a  true  contagion  for  their  demand. 

Most  of  the  seasonal  or  style  goods  change  annually  or  semi-annually. 

So  the  length  of  time  they  remain  in  market  is  much  smaller  compared  to 
that  of  new  product  lines.  Hence,  an  approximation  is  made  by  using  an 
(average)  constant  contagion  rate  throughout  the  season.  This  may  be  valid 
if  the  season  is  small  relative  to  the  variation  of  the  contagion  rate  or  if 
the  contagion  rate  is  changing  very  slowly. 

In  the  case  of  new  products  which  will  be  in  the  market  for  a 
considerable  length  of  time,  the  contagion  rate  is  approximated  as  a  step 
function  with  breaks  at  regular  intervals  of  time  identified  as  the  length 
of  the  review  period  T  .  Choice  of  T  is  made  institutionally  by  observing 
how  rapidly  the  contagion  rate  is  changing.  The  inventory  system  is 
reviewed  periodically  and  at  the  end  of  every  review  period,  a  decision  is 
made  about  the  best  stock  up  level  for  the  succeeding  period  depending  on  the 
present  inventory  level,  i.e.,  a  Periodic  Review  Policy  is  employed.  The 
unit  contagion  rate  is  not  assumed  to  be  the  same  from  one  period  to  another. 
Generally,  the  contagion  rate  will  be  less  in  each  succeeding  period  and  the 
successive  new  value  at  the  end  of  every  review  period  is  estimated  by  the 
knowledge  of  the  past  realizations  of  demand.  Under  this  assumption,  the 
unit  contagion  rate  at  time  t  ,  denoted  by  a(t)  ,  looks  as  follows  under  a 
particular  realization. 


FIGURE  1.1 


In  the  above  illustration,  the  product  has  three  review  periods  during  its 
transient  stage,  viz,  ,  T2  and  .  At  the  end  of  time 

+  T2  +  T^  ,  the  product  is  stabilized  in  the  market  and  the  contagion 
effect  vanishes. 

1.6  Assumptions  and  Notations 

Without  loss  of  generality,  the  beginning  of  each  period  is  taken  as 
time  zero.  For  the  present,  we  will  confine  ourselves  to  seasonal  goods  or 
to  the  first  review  period  of  new  product  lines,  where  at  time  zero  the 
demand  is  zero.^ 


This  is  not  true  for  succeeding  periods  of  new  product  lines  as  the  demands 
in  the  previous  periods  will  influence  the  demand  distribution  in  the 
succeeding  periods.  Though  this  does  not  complicate  our  analysis  too  much 
(luckily!),  we  postpone  its  discussion  to  Section  1.9. 
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Let  X  >  0  denote  the  initial  constant  demand  rate  component.  The 
unit  contagious  demand  rate  is  denoted  by  a  >  0  .  In  other  words,  a  is 
the  increase  in  demand  rate  for  each  past  demand.  N(t)  denotes  the  number 
of  demands  in  an  interval  [0,t]  of  length  t  .  The  probability  of  n 
demands  in  [0,t]  is  denoted  by  Pfl(t)  •  Let  T  denote  the  length  of  the 
review  period. 

Denoting  by  h  ,  the  length  of  a  very  small  interval,  we  make  the 
following  assumptions: 

Probabilities  of  positive  demand  occurrence  in  the  interval 
(t,t  +  h)  given  r  previous  demands  in  (0,t)  satisfy: 

(1.1)  P(N(t  +  h)  -  N(t)  ■  1/N(t)  *  r}  ■  (A  +  ar)h  +  o(h)  where  r  ■  0,1, 

and  o(h)  denotes  higher  order  terms  in  h  such  that  lim  »  0  . 

h+0 

(1.2)  P{N(t  +  h)  -  N(t)  >  2/N(t)  -  r}  -  o(h)  . 

(1.3)  P{N(t  +  h)  -  N(t)  -  0/N(t)  -  r}  -  1  -  (A  +  ar)h  +  o(h)  . 

1.7  The  Contagious  Probability  Distribution  for  Demand 


h 

I - 1 - 1 - 1 

0  t  t  +  h  T 

FIGURE  1.2 

In  the  above  illustration  of  the  review  period  [0,T]  ,  consider  an 
instant  of  time  t  c  [0,T]  .  Assume  the  interval  [t,t  +  li]  of  length  h 
satisfies  the  Properties  (1.1),  (1.2)  and  (1.3).  Hence,  we  can  write  the 
following  expression  for  the  probability  Pfl(t  +  h)  in  terms  of 
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probabilities  at  time  t  for  n  =  0,1,2 . 


For  n  *  0  : 


Po(t  +  h)  =  PQ(t)[l  -  Xh  +  o(h)]  . 


Taking  Po(t)  to  the  left  and  dividing  by  h  ,  this  reduces  to 


P  (t  +  h)  -  P  (t)  .  V 

- - 5  2  ■  -xp0(t)  *  It1  • 


Taking  the  limit  as  h  tends  to  0  ,  we  get 


P'(t)  -  -XP  (t)  . 
o  o 


Integrating  both  sides  and  using  the  initial  condition  Pq(0)  *=  1  ,  we  get 


PoU)  -  e 


-Xt 


Similarly,  for  all  n  >  1  : 

Pn(t  +  h)  »  Prob  [N(t)  ■  n  and  no  demand  in  (t,t  +  h)  ] 

+■  Prob  [N(t)  ■  n  -  1  and  one  demand  in  (t,t  +  h)] 

+  Prob  [N(t)  *  n  -  k  and  k(>  2)  demands  in  (t,t  +  h)]  . 


Writing  down  the  corresponding  probabilities  and  using  1.1,  1.2  and 
1.3,  we  get 

P  (t  +  h)  -  P  (t)[l  -  (X  +  ctn)h] 
n  n 

(1.4) 

+  P  ,  (t) [X  +  o(n  -  l)]h  +  o(h)  . 

n-i 

P^Ct)  to  the  left  and  dividing  by  h  ,  (1.4)  reduces  to 


Taking 
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p  (t  +  h)  -  p  (t)  0/hN 

— - - - 2 -  =  _(x  +  an)Pn(t)  +  [A  +  cx(n  -  DlP^U)  +  ^ 


Taking  the  limit  of  both  sides  as  h  tends  to  0  ,  we  get 


(1.5)  P'(t)  ■=  -(X  +  cm)P(t)  +  [X  +  a(n  -  1)]P  (t)  . 

n  n  n-i 


Define: 


p_L(t)  =  o 


for  all  t  >  0 


(1.6) 


P  (0)  -  1 
o 


P  (0)  =  0 
n 


for  all  n  >  0 


The  solution  of  the  above  differential  Equation  (1.5)  is 


P  (t)  ,  e"(X+na)t[x  +  a (n  -  1)]  f  e(X+na)tP  . (t)dt 
n  J  n-i 

(1.7) 

-(A+na) t 

+  ce 

The  general  solution  holds  for  all  values  of  n  *  1,2,  ...  and  the 
constant  of  integration  c  can  be  evaluated  using  the  initial  conditions 
given  in  (1.6).  It  has  already  been  derived  that 

P  (t)  *  Prob.  of  no  demand  in  (0,t) 
o 

-Xt 


Solving  (1.7)  for  n 
P:(t) 


1 

\ 

a 


,  we  get 


-Xt 

e 


+  ce 


- ( A+a ) 


t 


As  t  -*■  0  ,  P^ C t )  -*■  0  which  gives  c  =>  -X/a  .  Hence,  we  get  the 
probability  of  one  demand  in  (0,t)  as 
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,  x  X  -XtM  -at, 
PL(t)  -  -  e  [1  -  e  J  • 


n  =  2  and  n  =  3  ,  and  solving  for  c  using 


Similarly,  by  substituting  n  -  *  aim  u  » 

(1.6),  we  get  the  following  probability  expressions: 


P2(t)  -  ^ - 2 

2  2 1  a 


Mx  +  J*I  e"XtU  -  e~at] 2 


P3(t) 


\(\  +  a)(X  +  2a )_  e-Xt _  e"at]3  . 


Extending  the  results  to  n  ,  we  get 


Pn(t) 


•fa)  ...  [X  +  (n  -  l)al  e'Xt[i  -  e'at]n  • 


Taking  the  factor  a  out  from  each  term,  we  get 


(1.8)  Pn(t) 


-at.n 
1  -  e  ] 


Define  i-p  (a  constant) .  Since  X  and  a  ate  positive,  0  is 
also  positive!  But,  p  need  not  necessarily  be  an  integer.  Using  the 
substitution  X/o  -  0  ,  ve  can  rewrite  (1.8)  as 


+  „  -  IHp  1  n  -  2)  ■■■  (o  +  2  )(P  -X  (e"at)°  (1  -  e'at)n  . 


P  (t) 

n 


Using  the  well-kr  i  gamma  notation,  we  can  write 


P  m  «  r(p  +.n)  (e-at)P(1  .  e-otjn  . 

pn(t)  r(p)r(n  +  1) 


(1.9) 


r(p  +  n) 


To  simplify  the  writing,  we  will  henceforth  denote  r(p):(n  +  1) 


13 


(c  + " ' x)  •  »*•> 


r(n  +  1)  =  nr(n)  ;  T(4)  *  ^  . 


Thus,  the  contagious  distribution  becomes: 

(1.10)  PQ(t)  =  (p  +  "  1)(e“at)p(l  -  e“at)n  for  all  n  *0,1,2,... 


Equation  (1.10)  can  also  be  interpreted  in  some  sense  as  a  negative  binomial 
distribution  by  suppressing  time.  Then  Pn ( t )  will  be  the  probability  that 
exactly  n  failures  precede  the  pth  success,  if  p  is  an  integer  where 
prob.  of  success  is  given  by  p  *  e  at  >  0  and  prob.  of  failure 
~  q  *  1  -  e_at  >  0  for  all  t  >  0  ,  a  >  0  .  Hence,  Pn(t)  can  be  rewritten 
as 

(l.ii)  Pn(t)  c  +  n  "  1jppqn  for  a11  n  *  o.1.2.  ••• 

Note:  The  Equations  (1.10)  and  (1.11)  hold  for  n  =  0  also,  as  from 
(1.9), 


P  (t) 
o 


r(p  +  o)  o  p  -it 
r(p)r(o  +  i)  p  p 


Following  Feller  [A],  we  can  rewrite  the  binomial  coefficients  and  using  the 
fact  for  any  a  >  0 

(i)  -  <-»t + r x) • 


we  get 

(1.12)  Pn(t)  -  (’P)pD(-On 
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whore 

h)  c  (_Dn  — r.LP-±n)— 

\  n /  ^  T (p) F'(n  +  1)  ■ 

Tables  of  the  negative  binomial  distribution  are  available  in  [14],  Also, 
Taylor  [13]  shows  an  important  mathematical  equality  in  his  paper  that  the 
infinite  sum  of  negative  binomial  terms  can  bo  expressed  as  a  finite  sum  of 
(positive)  binomial  terms.  Hence,  the  latter's  table  can  be  used  to  find 
the  former. 

Feller  [4]  has  also  shown  that  Equations  (1.11)  and  (1.12)  represent 

an  honest  probability  distribution.  Hence,  for  a  fixed  "t,"  (1.10)  also 

represents  an  honest  probability  distribution. 

By  definition  the  generating  function  of  P^t)  »  denoted  by  Gfc(Z)  , 

00 

will  be,  G  (Z)  *  £  ZnP  (t)  .  It  can  be  easily  verified  that 

t  o  n 

n=0 


(1.13)  Gt(Z) 


where  p 


-at 

e 


-at 


1.8  Central  Moments  of  the  Contagious  Distribution 

Using  the  moment  generating  function,  the  mean  or  expectation  of  demand, 
for  fixed  t  ,  will  be,  by  definition 


m(t) 


A  6GtU) 
6Z 


Z-l 


Using  (1.13), 
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Hence, 

(1.14)  m(t)  ■  ^  ~~  f  )  . 

p  -at 

r  e 

Similarly,  taking  the  second  derivative  of  (1.13)  with  respect  to  Z  and 
taking  the  limit  as  Z  -*•  1  ,  we  get 

l  n(n  -  1)P  (t)  -  & -t V*  . 

„  n  l 

n=0  p 

By  definition,  the  variance  is  given  by 


V(t)  -  l  n2P  (t)  -  [m(t) ] 2  . 
n=0  n 

Hence, 

(1.15)  V(t)  -  &  =  "^.iat"  ~J'  • 

P  e 

1.9  Demand  Distribution  in  Second  and  Succeeding  Periods 

As  pointed  out  in  Section  1.6  in  the  case  of  new  product  lines,  the 
analysis  for  first  period  and  succeeding  periods  differ  mainly  because  of 
the  influence  of  first  period  demands  on  the  second  and  succeeding  periods, 
due  to  the  contagious  effect.  Luckily  enough,  the  analysis  is  not  com¬ 
plicated  too  much,  and  we  get  the  same  form  of  the  contagious  distribution 
as  derived  for  the  first  period,  except  for  a  change  in  parameter  depending 
on  the  number  of  demands  in  the  previous  periods.  The  reader  will 
immediately  notice  that  the  analysis  for  the  second  period  will  follow 
identically  to  third  and  succeeding  periods.  The  demand  distribution  for 
the  second  period  is  derived  assuming  the  number  of  demands  in  the  first 
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period  to  be  and  the  new  estimated  value  of  the  contagious  factor  be 

»2  •  Once  again  the  beginning  of  the  review  period  will  be  denoted  as  zero. 

Referring  to  Figure  1.2,  probability  of  a  demand  in  the  interval 
(t,t  +  h)  given  r  demands  in  (0,t)  will  be 

(1.16)  P{N(l  +  h)  -  N (t)  =  1/N (t)  =  r}  *  (X  +  +  c»2r)h  +  o(h)  . 

Since  we  know  a 2  and  ,  we  know  the  product  a2^l  a  "constant" 

demand  rate  and  can  be  added  the  other  constant  demand  rate  X  .  Let 

x2  ■  X  +  <*2N^  • 

Hence,  (1.16)  reduces  to 

(1.17)  P{N(t  +  h)  -  N (t)  *  1/N(t)  >=  r)  ■  (X  +  a2r)h  +  o(h) 
once  again 

(1.18)  P{N(t  +  h)  -  N(t)  >  2/N(t)  =  r}  =  o(h) 

and 

(1.19)  P{N(t  +  h)  -  N(t)  *  0/N(t)  »  r}  *  1  -  (X2  +  a2r)h  +  o(h)  . 

Immediately,  we  notice  that  the  simple  substitution  X2  =  X  +  N^oi2  ,  has 

reduced  the  Equations  (1.17),  (1.18)  and  (1.19)  identical  to  (1.1),  (1.2) 

and  (1.3).  Hence,  we  will  get  the  same  differential  Equation  (1.5)  except 

that  X  will  be  replaced  bv  X_  ,  which  will  lead  to  identical  solution  for 

X2 

P  (t)  .  Denoting  —  =  P-,  »  we  get 
n  a„  i 


(1.20) 


°2 

(2)  r(Po  +  n)  /  -a2C\  (  ”a2t\ 

’n  (t)  -  l  (pjr(n-T)-  U  /  U  -  6  I  11  --  ".I.*.  ■ 


Note  that  we  can  extend  these  results  to  any  succeeding  period.  For 

example,  let  denote  the  unit  contagion  rate  for  the  Nth  period.  Then 

N-l 

*  X  +  £  N^cx^  where  N^  is  the  number  of  demands  in  the  ith  period. 

x 

AN 

We  can  define  similarly  as  p^  *  —  ,  and  we  will  have 

N 


Pf>(0 


=n  +  °>  ( -v)°Y  -v)B 

r (n  +  l)  'e  /  V1  -  e  / 


for  all  n  »  0,1,2,  ... 


Thus,  we  notice  that  the  probability  distribution  for  demands  in  any 
review  period  depends  on  the  number  of  demands  that  occurred  in  previous 
periods.  It  is  interesting  to  note  how  the  central  moments  vary.  We  shall 
confine  ourselves  to  the  second  period  as  the  results  are  identical  for 


succeeding  periods. 


(t)  *  mean  number  of  demands  in  2nd  period 


,(■  • 


Since  P2  increases  linearly  as  the  number  of  demands  in  the  first  period 
(N^)  increases,  we  see  that  n^Ct)  t  ,  and  so  does,  1^(0  ,  the  variance, 


P2U  -  e 

V2(t)  “  -2a2t 

e 


1 . 10  Limiting  Distribution 


Since  it  is  postulated  that  the  contagion  rates  ...  are  strictly 
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decreasing,  it  will  be  important  to  know  the  limiting  form  of  the  contagious 
distribution  as  the  contagious  effect  vanishes.  It  turns  out  the  limiting 
distribution  is  Poisson. 

Theorem  1.1: 

The  limiting  distribution  of  P^(t)  »  given  by  Equation  (1.10),  tends 
to  a  Poisson  law,  as  a  ■*  0  . 


Proof : 

To  prove  that  Poisson  is  the  limiting  form  of  P^(t)  »  we  follow 
Feller's  [4]  approach,  using  the  generating  function.  From  (1.13),  the 
generating  function  of  P  (t)  is  given  by 


Gf  (Z) 


where 


P 


-at 

e 


p  +  q  =  1 


p  ■  X/a  . 


As  a-*0;p-»-'D,p-*-l  and  q  -*•  0  ;  let  pq  -*■  Xt  (fixed).  Taking  the 
limit  of  Gt(Z)  as  a  •»  0  ,  we  get 


lim  G  (Z) 
a-0 


lim 

Q-*KX> 


1  -  Xt/o  1° 

i  -  XZt/p 


lim  F(p) 

p-*x 


F(p) 


Xt/o  0 
XZt/p 


where 
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Taking  logarithms  of  both  sides,  we  get 

llm  loge  F(p)  -  lim  [l°«  <1  -  -  log  (1  -  \ZtM 

p-K»  p-XXt  L  '  ** 

Applying  L' Hospital's  Rule,  we  get 

lim  log  F(p)  »  -Xt  +  XZt  . 

p-foo 

Hence,  the  limiting  value  of  the  generating  function  as  a-*-0*eA,''‘'1' 
which  is  nothing  but  the  generating  function  of  a  Poisson  distribution. 
Hence,  as 


o  -  0  ,  Pn(t) 


Note: 

By  our  assumption,  in  a  multi-period  model,  the  contagion  rates  are 
strictly  decreasing  with  each  period  and  tends  to  zero  after  a  finite  number 
of  periods. 

Hence,  increases  with  N  and  tends  to  +»  in  a  finite  period. 

(N) 

From  the  proof  of  Theorem  1.1,  it  can  be  observed  that  P^  (t)  tends  to  a 
Poisson  distribution  after  a  finite  number  of  periods. 

1.11  The  Inventory  System 

In  Section  1.5,  it  was  shown  that  an  approximation  is  made  about  the 
contagion  rate  by  following  a  periodic  review  policy,  i.e.,  reviews  of  the 
inventory  system  at  stated  intervals  of  time  and  depending  on  the  inventory 
level  realized  at  the  beginning  of  each  interval,  the  best  ordering  policy 
is  chosen  for  the  succeeding  period.  Thus,  the  periodic  review  policy  leads 
to  a  multi-period  inventory  problem.  Because  of  the  computational  difficulty 
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of  not  having  an  explicit  expression  for  the  estimates  of  the  contagion 
factor,  the  problem  cannot  be  treated  as  a  multi-period  problem  but  only  as 
successive  single  period  problems.  At  the  beginning  of  each  interval,  using 
the  past  realization  of  demand,  the  new  value  of  the  contagion  factor  is 
estimated  which  gives  the  contagious  demand  distribution  for  that  interval. 

So  far,  nothing  has  been  assumed  about  the  stated  intervals  of  time. 
There  are  number  of  ways  to  determine  these  intervals  of  time.  One  way  is 
to  fix  this  review  period  institutionally  either  by  the  knowledge  of  past 
experience  or  arbitrarily.  In  this  case,  the  review  periods  may  be  of 
increasing  length  in  the  order  of  the  period,  to  facilitate  looking  at  the 
system  more  often  during  the  initial  periods,  which  one  might  normally  do 
since  the  "learning"  about  the  demand  behavior  is  limited  during  the  initial 
periods.  This  type  is  a  fixed  periodic  review  policy.  Chapter  2  essentially 
discusses  an  optimal  order  level  for  a  fixed  periodic  review  policy  for  a 
single  period  problem. 

Aga  in,  it  was  shown  in  Section  1.5  that  for  the  case  of  seasonal  or 
style  goods,  we  only  have  a  single  period  to  consider.  To  determine  this 
length  of  the  review  period  T  ,  one  may  follow  the  fixed  period-length 
policy  by  choosing  T  by  past  experience  or  arbitrarily.  Another  way  would 
be  to  seek  the  "best"  review  period  to  choose.  This  leads  to  an 
optimization  problem  with  two  variables,  namely,  an  order  level  and  a  review 
period.  This  type  is  an  optimal  period-length  policy,  which  can  be  solved 
as  a  "two-stage  optimization"  problem.  The  optimization  proceeds  by  first 
finding  the  optimal  order  policy  as  a  function  of  the  length  of  the  review 
period  and  then  selecting  an  optimal  review  period.  By  considering  the 
optimization  in  two  stages,  an  optimal  solution  for  a  fixed  period-length 

Chapter  4  discusses  the  estimation  procedures  for  the  parameters  of  the 
contagious  distribution. 


21 


policy  is  also  given.  This  is  discussed  in  Chapter  3.  This  is  the  first 
time,  to  the  author's  knowledge,  an  attempt  has  been  made  successfully  to 
find  an  optimal  review  period  for  an  inventory  policy  under  stochastic 
demand.  For  a  deterministic  demand,  we  have  the  classical  result  in  the 
Wilson-Harris+  lot  size  formula. 

Two  situations  arise  in  an  inventory  system  from  demands  which  occur 
when  the  inventory  is  zero: 

(i)  The  demand  is  backordered  and  supplied  at  the  start  of  the  next 
period,  i.e.,  backorders  are  filled. 

(ii)  The  demand  is  not  filled,  i.e.,  it  is  lost. 

Both  the  above  mentioned  situations  can  arise  for  seasonal  goods  as 
well  as  new  product  lines 

1,12  Costs  Associated  with  the  Inventory  System 

Three  types  of  costs  are  considered:  (i)  procurement  costs, 

(ii)  inventory  holding  cost,  and  (iii)  stock-out  costs. 

One  may  observe  in  inventory  literature  that  the  cost  function  is 
assumed  convex  without  specifically  giving  the  expression  for  it.  Instead, 
in  this  thesis,  each  of  the  costs  are  purposely  discussed  and  an  acceptable 
parametric  form  is  given.  In  practice,  this  facilitates  computing  these 
costs  more  easily  as  one  has  to  only  compute  the  various  parameters.  Also, 
this  gives  an  opportunity  to  make  a  parametric  study  of  the  optimal  solution 
and  see  how  sensitive  is  the  optimal  solution  with  respect  to  the  various 
cost  parameters. 


A  detailed  discussion  of  Wilson-Harris  lot  size  formula  can  be  found  in 
Hadley  [8] . 
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1.13  Procurement  Costs 

We  shall  begin  by  examining  the  procurements  costs,  which  can  be 

divided  into  two  parts.  First,  there  is  the  cost  of  the  goods  itself  which 

is  paid  to  source  of  production.  Then  there  is  the  costs  incurred  by  the 
inventory  system  itself  in  making  a  procurement.  These  may  include  costs 
associated  with  ordering,  bookkeeping,  transporting,  inspection,  testing  and 
so  on.  These  themselves  fall  into  two  classes,  one  which  depends  on  the 
amount  ordered  and  the  other  independent  of  the  amount  ordered.  Including 
the  cost  of  the  goods  itself  in  the  first  classification,  we  can  denote  the 
procurement  costs  as  the  sum  of  two  costs;  one,  which  depends  on  the  amount 
ordered,  denoted  by  c/unit  ordered,  and,  the  other,  independent  of  the  amount 
ordered,  denoted  by  k  ,  called  the  set-up  cost.  Hence,  the  total  cost  of 

placing  an  order  for  Q  units  will  then  be  k  +  cQ  . 

1.14  Inventory  Holdiug  Costs 

The  next  important  cost  to  consider  is  the  inventory  holding  cost. 
Included  in  these  are  the  real  out  of  pocket  costs  such  as  cost  of  insurance, 
taxes,  breakage  and  pilferage  at  storage  site,  warehouse  rental  and  costs  of 
operating  the  warehouse.  But,  the  most  important  cost  is  not  the  direct  out 
of  pocket  cost  but  an  "opportunity  cost"  which  is  the  cost  incurred  by 
having  the  capital  tied  up  in  inventory  rather  than  having  it  invested  else¬ 
where,  and  it  is  equal  to  the  largest  rate  of  return  which  the  system  could 
obtain  from  alternate  investments.  So,  we  shall  assume  that  the 
instantaneous  rate  at  which  inventory  carrying  charges  are  incurred  are 
proportional  to  the  investment  in  inventory  at  that  point  in  time.  Let,  I  , 
denote  the  holding  cost  in  $/unit  time/$  invested  in  inventory.  According 
to  Hadley  [8]  ,  the  reasonable  real  world  values  for  I  range  from  something 


like  0.15  to  0.35. 
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1.15  Stock-Out  Cost 

Back-order  costs  are  difficult  to  measure  since  they  can  include  such 
factors  as  loss  of  customers’  goodwill  (i.e.,  in  the  future,  he  may  take  his 
business  elsewhere).  Two  types  of  back  order  cost  for  our  inventory  system 
are  considered: 

(i)  A  back  order  or  shortage  cost  depending  only  on  the  amount 
backordered.  Denote  this  by  p/unit  backordered.  This  may  include 
the  cost  of  notifying  the  customer,  bookkeeping  of  the  amount  of 
back  orders  and  so  on. 

(ii)  A  variable  cost  depending  on  the  length  of  time  for  which  an 
order  remains  unfilled,  e.g.,  a  machine  shop  where  lack  of  parts 
keeps  the  machine  idle.  Denote  this  by  p/unit  short/unit  time  in 
back  order. 

1.16  Selection  of  an  Operating  Order  Policy 

The  purpose  of  this  chapter  is  to  find  an  "optimal"  ordering  policy  with 
the  help  of  a  mathematical  model  of  the  inventory  system.  So  far,  the  term 
"optimal"  has  been  used  loosely.  By  an  "optimal"  policy  we  mean,  the 
ordering  policy  that  maximizes  the  net  profit  or  minimizes  the  total  costs. 

In  some  situations,  like  the  ones  we  are  interested,  viz,  production  of 
consumer  products,  one  would  like  to  maximize  one's  net  profit.  In  some 
other  situations,  where  the  profits  are  always  negative  (e.g.,  post  office), 
one  may  like  to  minimize  the  cost.  In  some  cases,  both  these  considerations 
may  arrive  at  the  same  ordering  policy.  One  should  remember  that  the  net 
profit  or  total  costs  need  not  be  equivalent  to  a  strict  accounting  profit 
or  cost,  since  for  purposes  of  computing  optimal  ordering  policy  it  is  only 
necessary  to  include  those  costs  which  vary  with  the  operating  policy.  Costs 
which  are  independent  of  the  operating  doctrine,  like  the  cost  of  operating 
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the  information  processing  system  (which  includes  the  cost  of  making  an 

/ 

actual  inventory  count,  use  of  computer  or  the  cost  of  making  demand 
predictions),  need  not  be  included.  There  is  another  reason  why  the  profit 
or  cost  will  differ  from  what  would  be  computed  from  accounting  records. 
This  is  because  the  stock-out  costs  include  components  which  are  not  out  of 
pocket  costs,  like  the  opportunity  costs. 
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CHAPTER  2 

SINGLE  PERIOD  OPTIMAL  ORDER  POLICY 

2.1  Feed  Back  Control  Policy 

It  was  shown  in  Chapter  1  that  the  multi-period  inventory  problem  for 
the  case  of  new  product  lines  is  solved  by  considering  it  as  successive 
single  period  problems.  Thus,  this  chapter  will  be  devoted  exclusively  for 
finding  an  optimal  order  level  for  a  single  period  problem  under  a  fixed 
periodic  review  policy  where  the  demand  for  the  new  products  follow  a  conta¬ 
gious  law.  Because  of  the  fixed  periodic  review  policy,  the  review  periods 
are  chosen  by  the  knowledge  of  past  experience  or  arbitrarily. 

As  discussed  earlier,  the  inventory  model  calls  for  new  estimated 
values  of  the  parameters  of  the  demand  distribution  after  every  review  period. 
Hence,  it  is  assumed  that  the  values  of  the  constant  demand  rate  (X)  ,  the 
unit  contagion  rate  (a)  ,  the  number  of  demands  in  the  previous  periods 
(£  N^)  and  finally,  the  review  period  (T)  are  known.  This  will  completely 
specify  the  contagious  distribution,  given  by  Equation  (1.10). 

Since  a  single  period  problem  is  being  solved  at  the  end  of  any  review 
period,  the  past  experience  of  the  demands  in  the  previous  periods  is  used 
for  the  succeeding  periods.  Thus,  an  efficient  feed-back  control  policy  may 
be  determined  as  information  is  constantly  fed  back  into  the  inventory 
system  and  this  feed  back  is  used  in  determining  the  ordering  policy 
successively.  There  is  also  a  "learning  process"  associated  with  the 
operation  of  the  inventory  system. 

2.2  Notations 

We  shall  review  all  our  notations  and  symbols  used  so  far  and  introduce 


a  few  more  in  this  section. 


x  -  initial  inventory  level  before  reordering  goods. 

y  -  starting  inventory  level  after  reordering.  (Hence,  (y  -  x)  is  the 
amount  ordered.) 

T  -  review  period. 

N  -  amount,  demanded  during  the  review  period  [0,T]  . 

Note  that  N  is  an  integer  valued  random  variable  having  a  contagious 
probability  distribution 

Pn(T)  =  (°  +  ”  ^)ppqn  for  all  n  =  0,1,2,  ... 

where 

m  /p  +  n  -  l\ _ L(g...t.  n.) _ 

W  \  n  /  r(p)r(n  +  1) 

(ii)  p=»ea^,p>0,q>0  and  p  +  q  =  1  .  a  is  the  unit  contagion 
rate  for  the  current  period  under  consideration.  (a  >  0) 

(iii)  P^CT)  is  t^ie  probability  of  .  n  demands  in  [0,T]  . 

(iv)  p  =  X'/cl  and  A'  =  A  +  N'a  .  Where  A  is  the  constant  demard  rate 
for  the  current  period  and  N'  is  the  total  number  of  demands  in 
the  previous  periods.  Note  for  the  first  period  A '  =  A  since 

N'  =  0  . 

m(T)  -  mean  number  of  demands  in  the  current  period.  Note  m(T)  =  oq/p  . 
k -  set  up  cost.  (k  >  0) 
c  -  cost/unit  purchased.  (c  >  0) 

I  -  inventory  holding  cost  per  unit  time  per  $  invested  in  inventory. 

h  =  Ic  -  inventory  holding  cost  per  unit  tine  per  unit  held  in  inventory. 

(h  >  0) 

p  -  shortage  cost/unit  short.  (p  >  0) 
p  -  shortage  cost/unit  short/unit  time,  (p  >  0) 
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r  -  gross  revenue/unit  sold,  (r  >  0) 

2.3  Assumptions 
A-I: 

No  disposal  of  goods  is  allowed  at  the  end  of  the  period.  In  other 
words,  at  the  end  of  every  period  we  either  order  for  more  to  increase  the 
inventory  level  (=>  y  >  x)  or  order  nothing  and  stay  at  the  same  initial 
inventory  level  x(=>  y  *  x)  .  Hence,  y  >  x  always. 

A-II ; 

The  expected  net  revenue  from  unfilled  orders  are  not  included  in  the 
current  period.  Since  the  back  orders  are  supplied  only  at  the  beginning 
of  next  period,  the  revenue  from  filling  the  back  orders  is  included  in  the 
next  period  for  convenience.  This  assumption  is  not  correct  if  the  current 
period  is  the  last  period  since  there  is  no  succeeding  period.  This  will  be 
discussed  in  a  later  chapter.  (Refer  to  Section  3.8.) 

By  Assumption  A-I,  we  have  to  consider  only  two  cases: 

Case  (i)  *”0,y>Q,y>x 

Case  (ii)  x<0,y<0,y>x. 

The  relevant  costs  in  these  two  cases  are  discussed  using  Figures  2.1 
and  2.2  for  Case  (i)  and  Case  (ii)  respectively. 

2.4  Expected  Gross  Revenue  Function 

By  Assumption  A-II,  the  expected  revenue  corresponds  to  goods  sold  in 


the  current  period  only. 


Referring  to  Figure  2.1,  there  is  no  revenue  from  filling  last  period's 
back  order  if  x  is  nonnegative,  while  the  revenue  is  -rx  if  x  is 
negative.  Hence,  the  revenue  from  filling  back  orders  is  written  as 
-r  min  (x,0)  where  x  represents  the  algebraic  value  of  the  initial 
inventory  level  |x  =  oj  .  Similarly,  the  revenue  from  demands  in  the 
current  period  (Figure  2.1)  is  r  times  the  number  of  demands  for  realiza¬ 
tions  (a)  and  (b)  and  is  ry  for  (c).  Hence,  in  general,  the  revenue  from 

y  oo 

current  demand  is  r  l  nP  (T)  +  r  £  yP  (T)  .  Thus,  the  total  expected 

n-0  n  n*y+l  n 

y  ao 

gross  revenue  is  -r  min  (x,0)  +  r  l  nP  (T)  +  r  £  yP  (T)  .  By 

n=0  n  n-y+1  n 

changing  the  summation,  this  expression  becomes  -r  min  (x,0)  +  rm(T) 

00 

-  r  l  (n  -  y)P(T)  . 
n-y+1 

Case  (ii); 

Referring  to  Figure  2.2  and  Assumption  A-II,  the  gross  revenue  is 
r(y  -  x)  .  (Note:  x  <  0  and  y  <  0  .) 

2.5  Procurement  Cost  Function 

As  stated  in  Section  1.13,  this  includes  a  set-up  cost  and  the  cost  of 
purchasing  the  goods.  This  cost  is  the  same  for  both  Cases  (i)  and  (ii). 
Since  the  set-up  cost  is  positive  only  if  an  order  is  made,  a  delta  function 
is  introduced  to  take  care  of  the  case  when  no  order  is  made  (i.e., 
y  ■  x)  .  Hence,  the  set-up  cost  =*  k5(y  -  x)  where 

5 (y  -  x)  ■  1  if  y  >  x 


0  otherwise. 
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The  cost  of  goods  purchased  is  c(y  -  x)  .  Hence,  the  procurement  cost  is 
k6(y  -  x)  +  c(y  -  x)  . 

2.6  Expected  Inventory  Holding  Cost 
Case  (i): 

An  inventory  holding  cost  is  incurred  as  long  as  there  is  a  positive 
inventory  level  and  there  is  an  instantaneous  holding  cost  with  respect  to 
the  inventory  level  at  that  instant.  Considering  an  instant  "t"  e  [0,1]  , 
the  inventory  holding  cost  in  the  interval  [t,t  +  dt]  will  be 

f  (y  -  n)P  ( t ) d t  ,  and.  integrating  over  the  entire  interval  [0,1 ]  ,  the 
n=0  n 

expected  inventory  holding  cost  for  the  entire  period  is 

T 

P  (t)dt  . 
n 

Now,  the  integral  j  P^(t)dt 


(y  -  n) 


n=0 


J 


I 

f  ■  r.(c  +  nj.  .  (e-at)D(1  _ 

J  r(p) r (n  +  i)  >  u 


-at.n 
e  )  dt 


-at 


By  change  of  variables,  with  the  substitution  u  =  1  -  e  , 


*■  M 

/  V0dt  ■  / 


F(o  4-  n)  .0-1  n  du 

!-0)r(n  +  1)  (l  -  u)  u  T 


where 


q  -  1 


-aT 

e 
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The  integrand  is  an  Incomplete  Beta  Function  Ratio  (properties  of  incomplete 
Beta  functions  are  given  in  Appendix  II).  By  definition,  a  Beta  function 
with  parameters  (m,n)  is  given  by 


B(«,n) 


a. 

/m-1 
u  , 


U)n  ^du 


r  (m  +  n)  ’ 


and  the  incomplete  Beta  function  with  parameters  x  ,  m  and  n  is  given  by 


A 

Bx(mtn)  *  J*  u°  ^(1  -  u)n  ^du  for  x  <  1 


Define 


B  (m,n) 

ifco~  *  K 1 


for  x  <  1 


where  Ix(m,n)  denotes  the  incomplete  Beta  function  ratio.  Hence, 


(2.1) 


I_(n  +  l,p) 

P  (t)dt  -  -  -9-— - 

n  a  p  +  n 


and  the  expected  inventory  holding  cost  is 


Case  (ii) : 


I_(n  +  l,p) 

(y  - n)  — 


It  is  immediate  that  the  inventory  holding  cost  in  this  case  (Figure  2.2) 
is  zero,  as  there  is  only  a  negative  inventory  level  throughout  the  period. 


Pearson  [11]  has  tabulated  the  values  of  incomplete  Beta  functions  for 
various  parameters. 


1 
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2.7  Expected  Shortage  Costs 

As  stated  earlier,  there  are  two  kinds  of  shortage  costs,  one  which 
depends  on  the  amount  of  shortage  and  the  other  depending  on  the 
instantaneous  shortage  level.  Again,  Cases  (i)  and  (ii)  are  discussed 
separately. 


Case  (i):  (Figure  2.1) 

The  expected  shortage  cost  associated  with  parameter  "p"  is  nothing 
but  the  expected  shortage  level  at  the  end  of  the  period  times  the  cost 
per  unit  short  (p)  .  By  computing  the  instantaneous  shortage  cost  in  a 
fashion  similar  to  that  for  holding  cost,  the  total  expected  shortage  cost 
is 


oo  -  „  I_(n  +  l,p) 

P  l  (n  -  y)Pn(T)  +  £  I  (■>  -  y)  -VrS - 

n*y+l  n=y+l 


Case  (ii); 

Referring  to  Figure  (2.2),  the  total  expected  shortage  cost  is 


-py  +  pm(T) 


T 

-  P>T  +  P  / 

0 


m(t)dt 


which  reduces  to 


(p  +  p/a)m(T)  -  ppT  -  (p  +  pT)y  , 


using  the  fact 


m(t) 


p(l  -  e  at)  r  -at 
■  ~c  •  »le 


1]  ■ 
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2.8  Expected  Net  Revenue  Function 

The  expected  net  revenue  function  for  a  given  review  period  T  with 
parameter  y  is  denoted  by  Ti(y,T)  .  The  expressions  for  Tr(y,T)  for 
Cases  (i)  and  (ii)  given  separately  are: 

Case  (i): 


(2.2) 


y>0,x  =  0,y>x 

ir(y,T)  “  rm(T)  -  r  min  (x,0)  +  cx  -  k6(y  -  x) 

I_(n  +  l,p) 

+  —  f  (y  -  n)  — ^ — — - 

«  L  '  p  +  n 


-  cy 


n=0 


+  (P  +  r)  l  (n  -  y)P  (T) 
n=y+l 

„  I_(n  +  l,p) 

+  p/a  I  <n  -  y)  \  +  „  ' 

n“y+l 


P  (T)  and  I  can  be  computed  from  the  tables  available  in  [14]  and  [11] 
n  ~ 

q 

respectively. 

Case  (ii): 

y<0,x<0  ,y>x 
(2.3)  "(y,T)  «  -rx  +  cx  -  kS(y  -  x) 

-  [m(T)(p  +  p/a)  -  ppT  -  (p  +  pT  +  r  -  c)y]  . 

It  can  be  easily  verified  that  the  following  general  expression  for 

>  > 

both  Cases  (i)  and  (ii)  holds.  In  general,  >r  all  y=0,x-0,y>x 


(2.4) 


Tr(y,T)  *  rm(T)  -  r  min  (x,0)  +  cx 
-  [k«(y  -  x)  +  G(y,T) ] 
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where 


G(y,T)  =  cy  +  (p  +  r)  l  [n  -  min  (y,n)]P  (T) 

n=0  n 

„  I_(n  +  l,p) 

(2.5)  +  £  l  [max  (y,n)  -  n]  -a~  - 

O  A  p  »  II 


„  I_(n  +  l,p) 

+  p/a  l  [n  -  min  (y,n)]  -9p-  —  — - 


2.9  Properties  of  the  Cost  Function 

The  optimal  order  level  y  ,  that  maximizes  the  net  profit  Tr(y,T)  , 
is  obtained  by  minimizing  the  total  cost  function  k6(y  -  x)  +  C.(y,T) 
where  kf(y  -  x)  is  a  step  function  as  shown  below  in  Figure  2.3. 


FIGURE  2.3 

Hence,  the  main  interest  is  to  examine  the  properties  of  the  function 
G(y,T)  for  all  values  of  y  . 

Proposition  2.1; 

The  cost  function  G(y,T)  is  strictly  pointwise  convex  in 
y  e  {0,1,2,  >*}  ,  for  given  T  >  0  . 


Proof : 


Substituting  y  >  0  in  (2.5), 


(2.6) 


CD 

G(y,T)  =  cy  +  (p  +  r)  £  (n  -  y)P  (T) 

n=y+l 


(y  -  n) 


I_ (n  +  l,p) 

_i _ 

p  +  n 


oo 

+  p/a  l 

n=y+l 


(n  -  y) 


1  (n  +  1  ,p ) 

J _ 

p  +  n 


To  show  that  G(y,T)  is  strictly  pointwise  convex  in  y  ,  it  is 
sufficient  to  show  that  its  second  difference  (as  y  takes  integer  values 
only)  is  positive.  The  first  difference  of  G(y,T)  is 


AG(y.T)  £  G(y  +  l.T)  -  G(y,T) 


-(p  +  r  -  c)  +  (p  +  r) 


(2.7) 


n=0 


pn<T> 


I  (n  +  l,p) 

+  -  i  -*—z - 

a  Ln  p  +  n 
n*0 


K  I_(n  +  l,p) 

-  P/a  l  - 


n*y+l 


p  +  n 


The  second  difference  of  G(y,T)  is 


A2G(y,T)  *  AG(y  +  l.T)  -  AG(y,T)  -  (p  +  r)P  X(T) 


+ 


T_(y  +  2  ,p) 

f.h  i-pW _ 

\  a  /  p  +  y  +  1 


We  have  P  , .  (T)  >0  and  I  (y  +  2,p)  >  0  for  all  y  c  (0,1,2 . *} 

y  1  5 

and  T  >  0  .  Hence  A2G(y,T)  >  0  ,  implying  that  G(y,T)  is  strictly 
pointwise  convex  in  y  £  (0,1,2,  ...»  ™)  ,  for  given  T  >  C  . 
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Proposition  2.2: 

The  values  of  the  cost  function  G(y,T) [  and  G(y,T)|  (which 

y=0“  y=0 

represents  the  values  of  G(0,T)  reached  from  left  and  right  of  the  origin 
respectively)  are  equal. 

Proof : 

Putting  y  <  0  in  (2.5), 

(2.8)  G(y,T)  *  (p  +  r  +  p/a)m(T)  -  ppT  -  (p  +  pi  +  r  -  c)y  . 

Hence,  G(y,T)  is  linear  for  all  y  <  0  with  a  negative  slope  since 
r  >  c  .  Its  value  at  the  origin  is 

(2.9)  G(y,T)|  *  (p  +  r  +  p/a)m(T)  -  ppT  . 

y*0~ 

From  (2.6), 

(2.10)  G(y,T)  J  (p  +  r  +  p/u)m(T)  -  poT  . 

y=0 

By  (2.9)  and  (2.10),  the  values  of  G(y,T)  at  the  origin  from  the  left  and 
from  the  right  are  the  same. 

Proposition  2,3: 

G(y,T)  is  pointwise  convex  in  y  ,  for  given  T  >  0  . 

Proof : 

From  (2.8),  G(y,T)  is  linear  for  all  y  <  0  and  hence  pointwise 
convex  in  y  <  0  .  From  Proposition  2.1,  G(y,T)  is  strictly  pointwise 
convex  for  all  y  e  (0,1,2,  ....  °°}  .  Hence,  the  only  remaining  thing  to  be 


shown  is 
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4G(y,T)|  >  AG(y,T)| 

y=0+  y=0 


From  (2.8), 


AG(y  ,T) |  =  -(p  +  pT  +  r  -  c)  . 

y=0' 


From  (2.7),  putting  y  =  0  ,  we  get 


AG(y ,T) |  =  -(p  +  r  -  c)  +  (p  +  r)PQ(T) 

y=0 


I_(1,P) 

_i _ 


p 


I_(n  +  l,p) 

_i _ _ 

p  +  n 


using  (2.1)  that 


t 


AG(y,T)| 

y— 0 


(2.11) 


I_(n  +  l,p) 

P  (t)dt  -  -  -"Vt- -  .  we  have 

n  a  p  +  n 

-  -(p  +  r  -  c)  +  (p  +  r)PQ(T) 

T  a  T 

+  h  f  po(t)dt  -  ;  f  [  l  o  (  )  1 
o  o 

-  -(p  +  pT  +  r  -  c)  +  (p  +  r)Po(T) 


T 

+  (h  +  p)  /  PQ(t)dt  >  AG(y,T) 
0 


y*0 


Hence,  G(y,T)  is  pointwise  convex  throughout  in  y  for  all  T  >  0  . 

2.10  Critical  Order  Level 

Let  y  (T)  denote  the  critical  order  level  for  given  T  minimizing 
the  cost  function  G(y,T)  .  Since  G(y,T)  is  linear  with  negative  slope 
for  y  <  0  ,  it  is  immediate  that  the  global  minimum  of  G(y,T)  cannot  occur 
for  y  <  0  .  Since  G(y ,T)  is  strictly  pointwise  convex  for  y  in  non- 


negative  real  line,  yQ(T)  is  unique.  If  the  slope  (first  difference)  of 
G(y,T)  at  y  =  0+  is  positive,  then  the  minimum  of  G(y,T)  occurs  at 
y  *  0  .  If,  on  the  other  hand,  AG  (y  ,T)  at  y  =  0+  is  nonpositive,  then 
yo(T)  >  0  .  Hence,  it  is  necessary  to  consider  the  sign  of  AG(y,T)  at 

y  *  0+  .  Let 

a(T)  =  AG (y ,T) \  . 

y=0 


From  (2.11) 

T 

a (T)  -  -(p  +  r  +  PT  -  c)  +  (p  +  r)PQ(T)  +  (h  +  p)  J  PQ(t)dt  , 

0 

and  since  Po(t)  -  e  ,  it  follows  that 


Po(T)  . 

Thus,  If  a(T)  >  0  for  given  T  ,  then  yQ(T)  H  0  .  Otherwise,  y  (T)  >  0 
for  a(T)  ^  0  .  As  the  review  period  T  tends  to  zero,  a(T)  tends  to  a 
positive  value  c  .  Also,  a(T)  tends  asymptotically  to 
-(p  +  r  -  c)  +  —  -  pT  for  T  -*■  ,  which  is  negative  and  remains  so 

for  sufficiently  large  values  of  T  . 

Proposition  2.4: 

There  exists  a  unique  and  finite  T  >  0  such  that  a(T  )  *0  und 

o  o 

a(T)  0  for  T  <  T  . 

<  >o 

Proof : 


(2.12)  a(T)  •  -(p  +  r  +  pT  -  c)  + 


h  ■»  p 


+  p  +  r  - 


h-  ±  P 


From  (2.12) 
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a'(T) 


da(T) 

dT 


=  -  (P  +  r)]  AP„<T>  -  P 


and 


(T) 


-E 


h  +  p 

X 


i 


(p  +  r)|  X  Po(T) 


Consider  the  two  cases  where 


(a)  <  (P  +  r) 

A 

(b)  -  ;  -l  ‘  (p  +  r)  . 


Case  (a): 

Since  -  -  y  ^  <  (p  +  r)  ,  a'  (T)  is  strictly  negative  and  a(T)  is  a 
convex  decreasing  function  in  T  .  Hence,  there  exists  a  unique  Tq  such 
that 


a(T)  -  0  for  T  -  T  . 
<  >  o 


Case  (b): 

Since  —  ^  >  (p  +  r)  ,  a''(T)  is  negative.  Hence,  a(T)  is  strictly 

A 

concave.  Thus,  a' (T)  decreases  as  T  increases  and  becomes  negative 

after  some  T  and  stays  negative  thereafter.  Hence,  there  exists  a  unique 

and  finite  T  where  a(T  )  -  0  and  a(T)  >  0  for  T  <  T  . 

o  o  <  >o 

From  Proposition  2.4,  it  follows  immediately  that 


(2.13) 


y0(T)  i  0 

>  0 


for  all 
for  all 


0  <  T  <  T 
=  o 

ss  ~  «*\ 
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For  T  >  T 
=  o 


(2.14) 


the  critical  order  level  (integer)  y  (T)  is  sought  such  that 

AG(yo(T)  -  1,T)  <  0 
AG(yo(T),T)  >  0  . 


A  plot  of  the  function  G(y,T)  with  respect  to  y  ,  for  a  given  T  ,  is 
shown  in  Figure  2.4. 

2.11  Optimal  Order  Level 

Given  an  initial  order  level  x  ,  if  the  order  level  is  y  (T)  , 

(assuming  y  (T)  >  x) ,  then  the  total  cost  of  ordering  and  operating  at 
level  y  (T)  is  k  +  G(yo(T),T)  .  Instead,  if  the  order  level  remains  at 
x  ,  the  cost  is  G(x,T)  .  Hence,  the  optimal  policy  would  be  to  order  up  to 
yQ(T)  if  and  only  if,  G(x,T)  >  k  +  G(y  (T),T)  .  Otherwise,  no  additional 
stock  is  ordered.  Compute  a  level  s(T)  <  y^d)  ,  such  that 

(2.15)  G(s(T),T)  =  k+G(yo(T),T)  . 

Note  that  s(T)  may  be  positive  or  negative  for  given  T  . 

* 

Hence,  the  optimal  order  level,  denoted  by  y  (I)  ,  will  be 

Ix  if  x  >  s(T) 

. 

y  (T)  if  x  <  s(T)  . 

Denoting  y  (T)  by  S  ,  such  an  order  policy  is  known  as  an  (s  -  S)  policy, 
since  the  order  level  is  S  if  the  inventory  level  is  below  s  ;  otherwise, 
it  is  optimal  to  remain  at  the  level  x  . 

2.12  Algorithm  to  Find  the  Optimal  Order  Policy 


For  a  given  review  period  T  ,  the  optimal  order  level  is  found  from 
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y^(T)  by  using  (2.15)  and  (2.16).  Hence,  the  main  interest  is  to  have  an 
algorithm  to  find  yQ(T)  »  for  those  periods  T  which  are  greater  than  or 
equal  to  ,  since  yQ(T)  is  zero  for  all  values  of  T  which  are  less 

than  Tq  .  From  (2.14),  an  integer  y  (T)  is  sought  such  that 

AG(yo(T)  -  1 ,T)  <  0 
AC(yo(T),T)  >  0  . 

Assume  there  are  tv  >  nonnegative  levels  y^  and  y2(y2  >  y^  such  that 
AG(ylfT)  <  0  and  AG(y2 ,T)  >  0  . 

Note  that  is  the  critical  order  level  for  given  I  ,  if  y2  =  y^  +  1  . 

At  the  start  of  the  algorithm  the  interval  [ y ^ , y ^ ]  will  be  much  larger  than 
unity.  The  aim  is  to  reduce  the  interval  length  such  that  ultimately 
[y2  ~  y1  ]  equals  unity.  For  this  purpose,  the  well-known  bisection 
procedure  may  be  used.  Later  a  procedure  for  finding  the  initial  values  of 
y^  and  y^  will  oe  given. 

Bisection  Procedure: 

yl  +  y2 

(i)  Compute  the  median,  y  =  - — - . 

m  2 

(ii)  Calculate  AG(y  ,T)  using  (2.7). 

m 

(iii)  If  AG(y  ,T)  >  0  ,  then  set  y„  =  y  .  Otherwise,  set  y.  3  y 

m  2  m  i  m 

Now  we  have  a  new  and  shorter  interval  [ y ^  , y^  ]  • 

(iv)  If  the  length  of  the  interval  (y2  -  y^)  is  unity,  then  set 
y,,  =  yQ(T)  •  Otherwise,  go  back,  to  Step  (i). 

Computation  of  the  Initial  Values  of  y^  and  : 

Since  T  >  T  ,  we  know  from  Theorem  2.4,  that  AG(y,T)  <  0  . 

y=0 
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Hence,  we  set  =  0  .  Referring  to  Equation  (2.7),  it  is  possible  to 
find  an  M  ,  with  the  help  of  the  incomplete  Beta  function  tables  [11],  such 
that  I  (M  +  l,p)  -*•  0  .  In  which  case,  AG(M,T)  >  0  as  required.  For 

q 

example,  for  p  =  2  ,  I  _  -*■  0  as  M  -►  8  ;  I  0  ->  0  as  M  -*■  50  .  Hence,  we 
set  =  M  . 

2.13  Properties  of  the  Cricical  Order  Level 

It  will  be  of  interest  and  use  (refer  to  Chapter  3)  if  we  know  the 
properties  of  yQ(T)  as  the  review  period  varies  in  the  nonnegative  real 
line.  The  following  properties  of  y^fT)  with  respect  to  T  will  be  proved 
with  the  help  of  a  few  propositions: 


(1) 

y0<T) 

is 

nondecreasing  in  T  e 

[0,-)  . 

(ii) 

>o(r' 

is 

a  step  function  of  T 

e  [0,-) 

(iii) 

y0m 

is 

bounded  above  as  T  -* 

°°  . 

Define  a  function  FX,(T)  such  that 

N 


fn(i) 


(2.17) 


(p  +  r)  l  P  (T)  +  -  T  — “ — 7 - 

r  '  u  ^  Qt  w  o  +  n 

n>=0  n=0 


I  (n  +  l,o) 
00  —  ' 

p/»  I  J - 


n=N+l 


p  +  n 


for  all  N  -  0,1,2, 


Note  that  FX,(T)  equals  AG(N,T)  +  (p  +  r  -  c)  .  From  (2. 14),  N  is  the 
critical  order  level  if  and  only  if 


(2.18) 


FN_l(T)  <  (p  +  r  -  c) 
Fn(T)  >  (p  +  r  -  c)  . 


is  sought  such  that  F  (T) 


Hence,  for  a  given  T  ,  an  integer  N 


is  above 
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From  (2.1), 


T  I_(n  +  l,p)' 

—  —  —3 -  =  P  (T)  . 

dT  La  P  +  n  J  V  ' 


N 

l  [(p  +  n  -  DpaP^CT)  -  XPn(T)] 
n=0 

N 

+  (h  +  p)  l  P  (T)  -  p  for  all  N 

n=0  n 

To  treat  the  case  N  =  0  separately,  it  has  been  omitted  in  (2.20).  From 
(2.17),  putting  N  =  0  and  simplifying, 

A 

(2.21)  F  (T)  -  -pT  +  (P  +  r)P  (T)  +  [1  -  P  (T)  ] 

O  O  A  O 

(2.22)  F^(T)  =  -p  -  [X  (p  +  r)  -  (h  +  p))Po(T)  . 

Limiting  Values 

Next  is  to  examine  the  limiting  values  of  F^(T)  and  F^(T)  for  all 
N  =  0,1,2,  ...  as  T  0  and  T  -*■  =°  . 

(1)  T  -  0  : 

Fn(T)  -  (P  +  r)  for  all  N  =  0,1,2,  ... 

(2.23)  F'(T)  -  - 1 X (p  +  r)  -  h] 

o 

F’ (T)  -  h  for  X  >  1  . 

N  * 


Hence,  we  can  write 


F’(T) 


dFN(T) 

dT 


(2.20) 


(p  +  r) 


A  ft 
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(ii)  As  T  -»  °°  : 


F^(T)  -*•  for  all  finite  N  =  0,1,2,  .. 

and  for  p  >  0 

FjJ,(T)  -*■  -p  for  all  finite  N  =  0,1,2,  .. 

The  following  proposition  characterizes  the  behavior  of  F^(T)  and 
F^(T)  for  infinitely  large  N  . 

Proposition  2.6: 

The  critical  order  level  is  bounded  above  as  the  review  period  T  ->  00 

Proof: 

From  (2.20),  as  N  tends  to  , 

OO  ^ 

FJ(T)  -  (p  +  r)  l  l  (p  +  n  -  DpaP  X(T)  -  XPn(T)]  +  (h  +  p)  -  p 
n=0 

■  (p  +  r)pa  |p  +  ^  j  -  X(p+r)+h  =  h  for  all  T  >  0  . 

Hence,  for  infinitely  large  N  ,  F^,(T)  is  linear  with  positive  slope,  and 
from  (2.23)  FfjOO  >  (p  +  r  -  c)  for  all  values  of  T  >  0  and  for  N 
greater  than  some  bounded  integer.  Hence,  condition  (2.18)  cannot  be 
satisfied  by  an  infinitely  large  N  which  implies  that  y  (T)  is  finite 
for  all  values  of  T  . 

To  show  the  most  important  and  difficult  property  that  y  (T)  is  a 
nondecreasing  function  of  T  ,  it  has  to  be  proved  that  Fj,(T)  crosses  the 
"(p  +  r  -  c)  line"  at  most  once  for  all  finite  N  .  Through  the  following 


proposition,  a  much  stronger  result  is  proved. 
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Proposition  2.7: 


F^'  (T)  can  have  at  most  one  change  in  sign  in  the  range  T  e  [ 0 , °°) 


for  all  finite  N  *  0,1,2 . 


Proof: 


Consider  first  N  >  1  .  Rewrite  (2.20)  as 


F^(T)  =  (p  +  r)  l  [(p  +  n  -  lJpcxP^CT)] 
n=l 

N 

-  [X(p  +  r)  -  (h  +  p)]  £  P  (T)  -  p  for  all  N  >  1 

n=0  n 


Introducing  a  new  variable  u=q=l-ea  ,  write  F^(T)  as  a  function 

of  "u"  alone.  Note  as  T  varies  in  the  range  [0,“)  ,  u  varies  in  the 

half  open  interval  [0,1)  and  the  transformations  from  u  to  T  and  vice 

versa  are  one  to  one.  Hence,  P  (T)  reduces  to 

n 


P  (u)  =  c  (1  -  u)  u  where  c 
n  n  n 


T(o  +  n) 
rcP)r(n  +  i) 


and  F'(T)  is  transformed  to 
N 


F'(u)  *  ja(p  +  r)  l  (p  +  n  -  Dc^Cl  -  u)P+1un  1 
(  n=l 


N  ) 

-  [X(p  +  r)  -  (h  +  p)]  [  c  (1  -  u)Dun  -  p 


Taking  the  term  (1  -  u)D  outside  the  curly  brackets  and  noting  that 


(p  +  n  -  1)  c  =  nc  , 
n-1  n 
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^(u)  =  (1  -  u)P!a(p  +  r)p  £  cnutl  4  °(P  +  r)  l  nu  uD 
(  n=0  n=0  n 

N  A  N  \ 

-  a(p  +  r)  l  tic  -  [X(p  +  r)  -  (h  +  p)]  £  c  u"  - 

n=0  n  n=0  n  ) 


cancelling  the  like  terms,  this  simplifies  to 


(2.24) 


^j(u)  =  (1  -  u)D|(h  +  p)  £  cnun  -  (X  +  aN)  (p  +  r)cNuNj  -  p 


for  all  N  >  1  . 


For  a  check  note  that 


FjJj(O)  =  h  and  F^(l)  *  -p  . 

Since 

(2.25)  F”(T)  =  F^’(u)  ^  =  F”(u)a(l  -  u)  , 

a  >  0  and  (1  -  u)  >  0  for  all  u  e  [0,1)  ,  the  sign  of  F^' (T)  is  the 
same  as  that  of  F^' (u)  .  From  (2.24), 


F’’(u)  -  (1  -  u)1 


-  p(l  -  u) 


(h  +  p)  £  nc  un  ^  -  (X  +  aN) (p  +  r)Nc 


n»l 


n 


P-1 


N 


(h  +  p)  £  c  u11  -  (X  +  aN)(p  +  r)c 
A  n 


n*0 


,  N-l] 

,CNU  I 

n"| 

NU  I  - 


Taking  (1  -  u)P  ^  as  a  common  factor, 
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F”  (u) 


(1 


•  -it 


(h  +  p)  l  (p  +  n  -  l)c  .un  *(1  -  u) 
__i  n-i 


n=l 


-  [a  +  aN)(p  +  r)NcNuN  1(1  -  u)J 

-  jp(h  +  p)  l  cnu 

”] +  [c(i  +  aN)  (p  +  r)c^u^ 


(1  -  u) 


P- 


{(h  +  p)  l 
(  Ln=0 


N 


(p  +  n)c  un  -  l  (p  +  n  -  1) 


n=l 


-  (X  +  aN)  (p  +  r)Nc_,u^  "*■  +  (X  +  aN)  (p  +  r)Nck,u^ 

N  N 

.  N  ) 

-  p(h  +  p)  l  CQun  +  p(X  +  aN)  (p  +  r)c^u  j 

n=0  ) 


using  the  fact  that  (p  +  n  -  l)c  ,  =  nc  and  after  cancellation  of  like 
^  n-1  n 

terms, 


F”  (u)  *  (1  -  u)P  1|-(h  +  p)NcNuN  -  (X  +  aN)  (p  +  r)Nc^uN  1 

+  (X  +  aN)  (p  +  r)Nc^uN  -  p(h  +  p)cNuN  +  p(X  +  aN)  (p  +  r)cNuNj  . 


Hence, 


F' '  (u)  =»  c  (1  -  u)P  V  :{[(X  +  aN)  (p  +  r)  -  (h  +  p) ]  (N  +  p)u 

(2.26) 

-  (X  +  aN) (p  +  r)N}  for  all  N  >  1  . 

From  (2.26),  the  sign  of  F^'(u)  or  equivalently  the  sign  of  F^'(T) 
depends  only  on  that  of  the  terms  within  the  curly  brackets,  since 
(1  -  u)°  ^  ,  c^  ,  *  are  all  positive  for  u  s  [0,1)  .  The  sign  of  the 

terms  within  curly  brackets  in  (2.26),  depends  on  the  values  of  the 
parameters  and  there  are  only  two  possible  cases: 
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(X  +  aN)  (p  +  r)  >  (h  +  p)  . 


(True  in  practical  cases.)  In  this  case,  the  terms  within  the  curly 
brackets  can  have  at  most  one  change  in  sign  as  they  are  linear  in  u 
Hence,  F^'(T)  can  have  at  most  one  change  in  sign  in  the  range 
T  c  [0,-)  . 


(X  +  aN)  (p  +  r)  <  (h  +  p)  . 


In  this  case,  F',' (u)  <  0  for  all  u  £  (0,1)  =>  F'  '  (T)  <  0  for  all 
N  N 

T  e  (0,“0  .  Hence,  FjJ,'(T)  does  not  have  any  change  in  sign  in  the 


range  T  e  (0,“)  and  for  all  N  >  1  . 


Proposition  2.7  is  proved  for  all  N  >  1  .  For  the  case  when  N  =  0  , 


from  (2.22), 


(2.27) 


F”(T)  =  [X  (p  +  r)  -  (h  +  p)]XPo(T) 


This  implies  F' ' (T)  is  either  positive  or  nonpositive  depending  upon 

whether  X  (p  +  r)  -  (h  +  p)  >  0  or  X  (p  +  r)  -  (h  +  p)  <  0  respectively. 

Hence,  F' ' (T)  cannot  have  any  change  in  sign, 
o 

There  are  three  important  corollaries  which  follow  from  Proposition  2.7 
which  ultimately  prove  that  yQ(T)  is  a  nondecreasing  function  of  T  . 

Corollary  1: 

(a)  For  values  of  N  >  1  ,  F^Ci)  is  either  strictly  concave  throughout 
or  initially  strictly  concave  and  then  remains  strictly  convex  in 
T  in  the  range  [0,°°)  . 
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(b)  F  (T)  is  either  strictly  concave  throughout  or  strictly  convex 
o 

throughout  in  T  e  (0,®)  . 


Proof: 

Again  consider  the  two  cases  relative  to  the  magnitude  of  the 
parameters. 

(i) 


(X  +  aN)  (p  +  r)  >  (h  +  p)  . 

By  (2.27),  F"(T)  >  0  .  Hence,  F  (T)  is  strictly  convex  in  T  .  By 

o  o 

Proposition  2.7,  F^'(T)  can  have  at  most  one  change  in  sign  for  all 

N  >  1  .  From  (2.26),  note  that  F^' (T)  is  initially  negative.  Hence, 
F„(T)  is  either  strictly  concave  throughout  or  strictly  concave  at 
first  and  then  strictly  convex  throughout  as  T  ranges  [0,°°)  . 


(X  +  aN)(p  +  r)  <  (h  +  p)  . 

By  (2.27),  F”(T)  <  0  .  By  (2.26),  F ”  (T)  <  0  .  Hence,  F^ (T)  is 

strictly  concave  in  T  z  [0,®)  for  N  ■  0,1,2 . 

Corollary  2: 

For  all  values  of  N  z  {0,1,2,  ...}  ,  F^(T)  crosses  the 
"p  +  r  -  c  line"  at  most  once  and  from  above. 

Proof: 

Since  F^(0)  *  p  +  r  for  all  N  *  0,1,2,  ...  ,  F^(T)  starts  from  above 
the  "(p  +  r  -  c)  line."  Also,  F^CO)  *  h  for  all  N  3  1,2,  ...  and 
p^(co)  ■  -p  for  finite  N  >  0  .  Hence,  from  Proposition  2.7  and  Corollary  1, 


it  is  clear  that  once  F^(T)  becomes  negative,  it  stays  negative  thereafter. 
Hence,  Corollary  2  is  true.  For  those  N  which  are  greater  than  some 
bounded  integer,  F  (T)  increases  continuously  and  hence  will  never  cross 
the  "p  4  r  -  c  line." 


Corollary  3: 


The  solutions  Tq,T^, 


T  (where  N  is  the  limiting  value  of 

N-l 


y  (T))  obtained  by  equating  F  (T)  =  p  +  r  -  c  for  N  =  0,1,2,  ...,  (N-l) 

O  N 

respectively  ire  unique  and  are  strictly  increasing. 


Proof : 


Using  (2.18)  and  Corollary  2  and  the  fact  that  N  is  the  limiting  value 

of  y  (T)  ,  it  is  immediate  that  F^(T)  crosses  the  "(p  +  r  -  c)  line" 

only  once  for  \  -  0,1,2 . (N-l)  .  Hence,  the  solutions 

Tn,T. ,  ...,  T  are  unique.  By  Proposition  2.5,  it  is  immediate  that 
U  '  N-l 

T0  <  T1  <  T2  <  <  T-  ,  * 


Example : 


Figure  2.5  illustrates  the  shape  of  F^(T)  with  respect  to  T  e  [0,®)  , 
for  the  case  when  (\  +  ctN)  (p  +  r)  >  (h  +  p)  ,  for  values  of 

N  -  0,1,2,  . . . ,  6  . 


Proposition  2.8: 

The  critical  order  Level  is  a  nondecreasing  step  function  of  the  review 
period,  continuers  from  the  right  with  finite  saltus  at  points 


T 


0 


T 

N-l 


Proof : 


From  Corollary  2  and  Proposition  2.5,  it  follows  for  a  given  T  =  T  , 


FIGURE  2.5 


that  if  the  critical  order  level  In  y^tT)  ,  then  for  all  T  >  T  , 

condition  (2.1H)  can  only  be  ant Uf led  by  an  N  y  (T)  ,  Hence, 

*  o 

•  • 

y0(T)  >  ye(T)  for  all  T  >  T  .  The  rest  of  the  proposition  follows  from 
Corollary  3  and  Equation  (2. IB).  Titus, 

yQ(T)  0  for  0  <  T  <  Tq 

y0(T)  -1  lot  T0  <  T  <  T, 


y  <T)  •  8  -  1  lor  T  <  I  <  T 
°  N-2  N-l 

y  (T)  -  N  for  T  >  T 

0  "  N-l 


Example : 


Figure  (2.6)  illustrates  the  plot  of  yQ(T)  against  T  for  the 
example  given  in  Figure  2.S. 


FIGURE  2.6 
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CHAPTER  3 

S1NCLK  PERIOD  MODEL  FOR  SEASONAL  OR  STYLE  GOODS 

3.0  Introduction 

As  discussed  In  Section  1.11,  for  the  esse  of  sessonnl  or  style  goods, 
either  s  fixed  period-length  policy  or  sn  optimal  period-length  policy 
may  be  followed.  In  the  former  esse,  sn  optimal  order  policy  Is  sought  by 
fixing  the  review  period  Institutionally,  while  In  the  latter  case  sn 
optimal  order  policy  as  well  as  an  optimal  review  period  are  sought  which 
leads  to  an  optimization  problem  with  two  variables.  The  optimization  Is 
carried  out  In  two  stages  by  first  finding  the  optimal  ordering  policy  as  a 
function  of  the  length  of  the  review  period  and  then  selecting  an  optimal 
review  period.  Thus,  an  optimal  solution  for  both  fixed  and  optimal  period- 
length  policies  are  given.  Only  a  single  period  problem  is  solved  for 
seasonal  goods  since  in  general  the  season  is  not  longer  than  a  year. 

Here  again  two  situations  may  arise  from  demands  which  occur  when  the 
Inventory  is  zero: 

(1)  The  demand  is  not  filled,  i.e.,  it  is  lost.  The  cost  of  losing 
these  demands  may  be  estimated  by  the  stock-up  cost  parameters 

A 

p(>  0)  and  p(>  0)  as  explained  in  Section  1.15.  If  the  cost 
of  "lost  sales"  is  ignored,  then  both  p  and  p  will  be  zero. 

(ii)  The  demand  is  back  ordered  and  supplied  at  the  end  of  the  season. 
In  general,  it  is  not  practical  to  follow  this  policy  for 
seasonal  goods. 

3.1  Optimal  Order  Level  for  a  Given  Review  Feriod 

Confining  the  analysis  to  the  case  when  "lost  sales"  policy  is  followed, 
it  is  immediate  that  all  the  results  of  Chapter  2  h:>ld  by  making  the  initial 


s  t 

order  level  (x)  noro,  Thu*  the  critical  order  level  yQ(T)  satisfies, 

AC(ye(T)  -  1,T)  <  0 
AC(y0(T),T)  >  0 

where  AG(y,T)  la  given  by  (2.7).  The  optimum  order  level  y*(T)  la, 

yo(T)  If  k  +  C(ye(T),T)  <  C(0,T) 

*  . 

y  (-)  ■< 

,  0  otherwise. 

* 

Since  a  alngle  period  problem  la  solved  for  seasonal  goods,  y  (T)  equals 
zero  implies  no  business  is  done.  Hence,  to  be  in  business,  a  positive 
Inventory  level  must  be  maintained  at  the  beginning,  tor  this  to  be  true, 
the  condition  k  +  G(yo(T),T)  <  G(0,T)  must  be  satisfied.  Proposition  3.7 
(refer  to  Section  3.5)  proves  that  there  exists  a  lower  bound 
T^(>  Tq  and  finite)  on  the  review  period  T  ,  above  which  y  (T)  is 
yo(T)  (>  1)  .  Thus,  only  those  values  of  T  which  are  greater  than 
will  be  considered  for  finding  an  optimal  review  period. 

3.2  Order  Level  Optimized  Net  Revenue  Function 

The  order  level  optimized  net  revenue  function  is  defined  to  be  (see 
(2.4)) 

tt(T)  -  7r(y  (T),T) 

(3.1) 

-  rrn(T)  -  k  -  G(yo(T),T)  for  all  T  e  (T0,~) 


and  (see  (2.5)) 
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C(y  <T),T)  •  cy  (T)  ♦  (p  4  r)  J  (n  -  y_(T) 1 P  (T) 
°  #  n-v  >)+l  °  " 


y.(T) 


l_(n  ♦  l,o) 


0.2) 


* :  j0  !».«)  -  -Vt  n 


1  (n  +  1,0) 


5  * ... ,n  -  yo<I)l  _1~ 


n»yo(T)+l 


From  (3.1)  and  (3.2),  It  la  observed  that  the  order  level  optimized  net 

revenue  ia  a  function  of  the  review  period  alone.  Hence,  an  optimal  review 

period  can  be  sought  by  maximizing  *(T)  with  respect  to  T  .  Recall  the 

fact  (tropoaitlon  2.8)  that  yQ(T)  is  s  step  function  of  T  ,  having  jumps 

at  points  Tn,T1#  ...,  T  where  at  time  T„  ,  y  (T)  Jumps  from  a  value 

0  1  N-l  N  0 

N  to  a  vlue  N  +  1  ,  being  continuous  from  the  right  and  N  is  the 

limiting  value  of  yQ(T)  . 

Since  the  expressions  for  r(T)  involves  yQ(T)  ,  it  will  be 

interesting  to  see  whether  tt(T)  has  also  jumps  at  T0,T.,  ...,  T  .  It 

N-l 

turns  out  that  *(T)  is  continuous  for  all  values  of  T  . 


Proposition  3.1: 

The  order  level  optimized  net  revenue  function  is  a  continuous  function 
of  the  review  period. 

Proof; 

From  Proposition  2.8,  it  follows  that  y^CT)  is  constant  during  the 
intervals  (T  ,T  ),[T  ,T  ),  ...,  T  ,T  T  .  Denote  by  r,  (T)  , 

U  L  N-2  N-l/  L  N-l  / 

the  value  of  the  function  it(T)  during  the  interval  [T  .  ,T.,)  where  y  (T) 

N-l  N  O 

is  equal  to  N  .  With  the  notation  T  *  ®  ,  N  varies  from  1,2,  ...,  N  . 

N 


Hence , 


5* 


0.3) 


tM<T)  •  r«(T)  -  k  -  C(N,T) 


for  all  Tw_j  <  T  < 
and  N  •  1,2,  . . . ,  R 


vhara 


C(N,T)  -  cN  +  (p  +  r)  l  (n  -  N)P  (T) 

n-N+1  n 


0.4) 


I.(n  +  l,p) 


+  -  l  (N  -  n)  — - — 

*  o=0  0  +  n 


+  5  !  <n 

a  n-N+1 


I  (n  +  l.o) 

N)  -i-r - 

p  +  n 


Ic  la  clear  from  (3.3)  and  (3.4)  that  itm(T)  la  continuous  In  the  Interval 

N 


<  T  <  TN  for  all  N  -  1,2,  ...,  N  .  The  only  thing  remains  is  to 
that  the  jump  at  point  is  zero.  Denote 


show 


4W  ■  WV  -  VV  • 


By  (3.3),  the  Jump  at  T  •  ,  is 


4»h<tn)  -  _Gh+1  +  VV 

■  -4vv  ■ 


In  Section  2.13  (Corollary  3),  it  was  shown  that  satisfies  AG^T^)  -  0 
Hence , 


AW  E  0  for  a11  N " 1*2*  •••»  <N_1^  • 


3.3  Limiting  Values  of  tt (T) 


The  next  question  would  be  to  examine  the  limiting  values  of  ir(T)  as 


T  tends  to  zero  and  Infinity. 
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T  -  0  J 

When  the  review  period  tends  to  zero,  yQ(T)  ■*  0  by  (2.13) 


<1 

J  un(l  -  u>°  *du 
I_(n  ♦  l,p)  »  - - 

£  Un(l  -  u)°  *du 
0 


where 


q  -  1  -  e 


-aT 


As  T  0  ,  q  -*•  0  .  Hence,  I_(n  +  l,p)  0  .  Also,  m(T) 

\ 

Pn(T)  0  for  all  n  >  1  .  Hence,  from  (3.1), 


11m  ir(T)  ■  0  . 
T-0 


)  T  ■»  °°  : 

Now,  as  T  *  ,  yQ(T)  -►  N  which  is  the  finite  upper  bound 
Also,  as  T  ->■  ®  ,  q  -*■  1  .  Hence, 

I_(n  +  l,p)  •*  1  . 

q 

Also,  P  (T)  ■*  0  for  all  n  >  0  and  m(T)  00  .  Now, 

*  n  • 

11a  ir(T)  -  lim  [rm(T)  -  k  -  G(yo(T),T)]  . 

X-**® 


Now, 


0  and 


for  y0(T)  . 


From  (3.2), 
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G<y  (T),T)  -  cy  (T)  +  (p  +  r)  f  In  -  yft(T)JPn(T) 
°  0  n-yjT)+l  0  n 


h  y°!T)  Vn  +  1't>) 


0  '0 

O-0 


P  ♦  n 

I-(n  +  l(p) 


♦  *  I  [n  -  y-(T)l  —  — - 
*  n-y0(T)tl  0  c  +  " 


Changing  the  order  of  summation  of  the  terms  containing  (p  +  r)  ,  we 
get. 


C(yo(T),T)  -  cyo(T)  +  (p  +  r)  £  [n  -  yo(T)]Pn(T) 

n-0 


y0(T) 

-  (P  +  O  l  [n  -  y0(T)]Pn(T) 
n-0 


y0(T) 


•  w  I-(n  +  1,0) 

+  *  l  p  +  n  — 

n-0 


*  ”  I-(n  +  l,p) 

J  l  [n  -  y  (T)J  q  ■  - 

°  n-y0(T)+l  °  P  +  n 


Using  the  fact  that  T  P  (T)  ■  1  and  l  nP  (T)  -  m(T)  ,  we  get, 

“  n  _  n 

n-0  n-0 


G(yQ(T),T)  -  cyo(T)  +  (p  +  r)[m(T)  -  yQ CT) ] 


y0(T) 

+  (p  +  r)  l  [y„(T)  -  n]P„(T) 
n-0 


(3.5) 


.  VT)  I-(n  +  l,p) 

+  !  l  lyon)  - ni 
n-0 


+  £  I  t"  -  y„tt» 


I-(n  +  l,p) 

_S - 


p  +  n 
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Pro®  (3.5),  «•  T  •*  *  , 


cy0<T)  ♦  cN 

(p  +  r) (m(T)  -  y0(T))  -  • 


h 

a 


ya(T) 

l 

n-0 


Iyo(T)  -  nj 


I_(n  +  l,o) 

J _ 

0  +  n 


IN  -  n) 


1 

o  +  n 


<n 


and  positive 


*  m 

£  l 

“  n-y0(T)+l 


In  -  yo(T)J 


I_(n  +  l,p) 

_g _ 

0  +  n 


£ 

a 


l  In  -  N) 


p  +  n 


n*N+l 


Hence,  11m  G(y  (T),T)J  -  •  .  Thus,  lim  rr (T)  •  »  -  »  which  Is  an 
T-h»  °  T-*<» 

Indeterminate  form.  L' Hospital's  Rule  may  be  used  to  find  the  limit.  Since 
the  limiting  value  of  ir(T)  is  not  in  the  standard  form  »/«  or  0/0  for 
direct  application  of  L'Hospital's  Rule,  transform  the  expression  in  such  a 
way  so  as  to  get  “V00  form. 

■*  _aT 

Multiply  the  expression  (3.1)  throughout  by  p  ■  e  to  obtain 

ir(T)  “  7  [rpq  -  pk  -  pG(yo(T),T)]  . 

P 


This  expression  can  be  rewritten,  using  (3.5),  as 


"(T)  -  7  \-ppq  -  kp 
P 


+  P 


-  P 


y  (T) 


(p  +  r  -  c)y  (T)  -  (p  +  r)  l  [y  (T) 

n-0 


yo(T)  I_(n  +  l,p) 

h  r  r __  />r>\  _i  SL 


n]Pn(T) 


n*0 


p  +  n 
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l  I»  - 


I_(n  +  l,p) 

__a _ _ 


n=yo(T)+l 


>.»!  p  +  „ 


1/P 


r  (Qx  +  q2  +  Q3  +  q4) 


i 

j 


where 


*1 


q2  *  p 


q3  *  -p 


-  kp 

y0(T) 

" 

+  r  - 

c)y0<T) 

-  (P  +  r)  I 

n=0 

[yo 

(T)  -  n]Pn(T) 

y0(I) 

I_(n  +  l,p) 

l 

ly„(T)  - 

n]  — - - 

p  +  n 

n=0 

00 

I_(n  +  1, 

P)1 

1 

l  [n  -  y  m]  -S— n 

/  rr.  \  i  i  O  p  *t  Tl 


n=yo(T)+l 


1/p 

I  (n  +  l,p) 


As  T  -*•  00  :  y  (T)  ■+  N  ;  — — — 

Jo  ’  p  +  n 


'  *  TT7  :  Qx  -  -Vt  ;  Q2  *  o 

and  -►  °°/®  .  So  apply  L'Hospital's  Rule  for  .  Hence, 


lim  [Qj 
T-»« 


lim 

T— 


-P  I  In  -  y  (T)]P  (T) 

(1  ;-rl  0  " 

_ _Q _ 

a/p 


Changing  Che  order  of  summation. 


lim  [Q,] 

X-*«> 


lim 


-pm(T)  +  pyQ(T) 


.  VT> 

P  l  [y  <T>  -  n]Pn(T) 
n*=0 


a/p 
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*-  -  .  VT> 

.  ££a  .  £  P  l  (y  (T) 

0  *  n-0  0 


n)Pn(T)  -  py0(T) 


•  -pp/a  . 

Hence,  the  numerator  of  *(T)  ,  which  ie  nothing  but  (Q^  +  +  Q^) 

A 

tends  to  a  constant,  -(p  +  p/a)p  ,  while  the  denominator  p  tends  to  aero. 

Hence,  lim  *(T)  •  . 

T-» 

3.4  Derivative  of  it(T) 


In  expressions  (3.1)  and  (3.2),  n(T)  Involves  the  step  function  yQ(T) 

which  does  not  have  left-hand  derivatives  at  points  T.,T  ,  ...,  T 

0  1  N-l 

Hence,  fixing  the  value  of  yQ(T)  ■  N  where  it  is  constant  for  all 
Tn1  <  T  <  TN  ,  its  derivative  may  be  examined  with  respect  to 
T  e  [T^.y  .  From  (3.3)  and  (3.4), 


nN(T)  -  rm(T)  -  k 


-  cN  -  (p  +  r)  l  (n  -  N)P  (T) 
n-N+1  n 


I  (n  +  l,p) 


I  (n  +  l,p) 


-  -  l  (N  -  n)  -  a-~ T  - J  l  (n  -  N)  . 

a  n-0  P  +  n  0  n-N+1  p  +  ° 


for  all  T„  ,  <  T  <  T„  and  for  all  N  -  1,2,  ...,  N  . 

n-l  -  N 


(Note:  T  -■+*>.) 
N 


*»«> 


Using  (2.1)  and  (2.19), 


dT 
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<(T>  -  ~  -  (p  +  r)  J  (n  -  N){(p  +  n  -  i)pa?  (T)  -  XP  (T)] 

"  p  n-N+1  n_1  n 

N  .  - 

-  h  l  (N  -  n)P  (T)  -  p  J  (n  -  N)P  (T) 
n-0  n  n-N+1  n 

change  Che  order  of  summation  for  p  term  to  obtain, 

C(T)  -  “  -  (h  +  p)  l  (N  -  n)P  (T)  +  p[N  -  m(T)] 
p  n-0 

Ot' 

<3,6)  -  (p  +  r)  l  (n  -  N)  f  (p  +  n  -  l)pcPn  (T)  -  AP  (T)] 

n-N+1  n_i  n 

for  all  <  T  <  TN  And  for  all  N  -  1,2,  ...,  N  . 

Proposition  3.2; 

For  all  N  ■  1,2,  ...,  N  ,  the  function  ir^(T)  is  continuous  in 
T  c  lTN_i*TN)  *nd  has  a  positive  Jump  at  the  boundary  where  yo(T) 

changes  its  value  from  N  to  (N  +  1)  . 


Proof : 

From  (3.6),  it  is  obvious  that  ir^(T)  is  continuous  in  the  interior. 

The  only  thing  is  to  prove  the  jump  at  T  ■  TM  is  positive. 

n  v" 


AirN(V 


-  jump  of  the  derivative 


*;(T) 


at  T  - 


TN< 


N+l  N 


<T„)  -  *'( T„) 


N  N 


From  (3.6),  the  above  will  be  equal  to 
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Al,N(V  "  P  "  <h  +  P>  I  (N  +  1  -  n)P 


"N+l 

l 

_  n=0 


W  -  I  w  -  ">pn<yl 

n-0  J 


) 

_n 


-  (p  +  r)  l  (n  -  N  -  l)](p  +  n  -  D^oP^  (TN)  -  AP^)) 


n=N+l 


=N+2 


(n  -  N) [  (p  +  n  -  l)pcPn_1(TN)  -  AP^)] 


Expanding  the  terms  under  summation  and  after  cancellation, 


4W  ■  p  -  <»  +  p)  l  W 

n=0 


N 


-  (p  +  r)  l  [<P  +  n  -  DpaP^Cy  -  AP^)] 
n*0 


The  above  expression  and  (2.20)  are  identical  except  for  a  sign  change. 

Hence,  A7r'(TXT)  *  -F'(T„)  .  From  Corollaries  2  and  3,  in  Chapter  2, 

N  N  N  N 

(3,7)  FjJ(T)  <  0  for  all  T  > 


since  F.,(T)  crosses  the  "(p  +  r  -  c)  line"  at  TXI  and  continues  to 
(strictly)  decrease.  Hence,  An^(T^)  >  0  . 

It  will  be  easy  to  investigate  the  behavior  of  it^(T)  if  there  are 
finite  limits  in  the  summation.  Hence,  by  changing  the  infinite  limits  in 
the  terms  corresponding  to  "(p  +  r)"  and  "p"  in  expression  (3.6)  to 
finite  limits,  (3.6)  can  be  rewritten  as. 


tt'(T)  -  ~  -  (h  +  p)  l  [N  -  n]P  (T)  +  p [N  -  m(T)] 


N 


n=0 


+  A (p  +  r) 


l2*  -  l  nP  (T))  -  n(i  -  l  P  (T)j 
\  p  n*0  /  '  n=0  > 


-  a(p  +  r)p 


l  (n  -  N)[pPn_1(T)  +  (n  -  DP^T)] 


n-N+1 
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L(T)  =  -a(p  +  r)p  l  (n  -  N)[pP  ,<T)  +  (n  -  1)P  ,(T)] 

|_n=N+l 

-  a(p  +  Op|n|p  -  P  l  Pn_i(T)j  +  N^m(T)  -  ^  (n  -  DP^O^ 

+  r)p  £ 

+  [(n  -  l)2  +  (n  -  DJP^id)) 


'] 


-  o(p  +  r)o|  l  jp[(n  -  DP^iCT)  +  Pn-1(T)] 
L  n-N+1  * 


Once  again  changing  the  infinite  limits  in  the  summation  to  finite  limits 
and  using  the  fact  that 


l  n2Pn(T)  -  V (T)  +  [m(T)]2  - 
n-0  P 


we  have 


(  N 

L(T)  -  a(p  +  r)p  N[p  +  m(T)]  -  N  l  (p  +  n  -  l)Pn_1(T) 
l  n-1 


N  N 

■  p  m(T)  +  P/p  -  P  l  nP  ,(T)  -  l  n(n  -  1)P  . 

n-1  n  1  n-1 

— 


(T) 


+  [m(T)]2  {  , 

P  J  ) 


1  •  ©  •  t 


(  N 

L(T)  -  -a(p  +  r)p  J  — N t p  +  m(T)]  +  l  [N  -  n]  (p  +  n  -  l)Pn_^(T) 
(  n-1 


4fP  +r  +  4  +  rl! 

p  L  q  P  P  J  ) 


Hence , 
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in  turn,  leads  to  the  investigation  of  irf(T)  .  Since,  the  derivative  is 
defined  only  in  the  intervals  where  yQ(T)  is  constant,  each  of  the  intervals 
is  investigated  separately.  Note  that  there  are  only  a  finite  number  of 
intervals.  As  a  matter  of  fact,  there  are  only  N  intervals  since  yQ(T) 
is  bounded  above  by  N  . 

For  all  fixed  yQ(T)  *  N  in  the  range  TN_^  <  T  <  TN  ,  expression 
(3.8)  may  be  written  as, 


(3.10)  and  (3.11)  hold  for  all  N  ■  1,2,  ...,  N  .  Now  investigate  S^CT) 
with  respect  to  T  e  [TN_^,TN)  •  Once  again,  make  use  of  the  transformation 
u  -  5  -  1  -  e-aT  and  Pn<T)  -  cn<l  -  uj^u"  where  ■  r{$(*  fjj  • 

Since  u^  *  1  -  e  ,  the  problem  reduces  to  examining  S^(u)  with  respect 
to  u  e  [u^J_1,uN)  .  Thus,  from  (3.10), 

SN(u)  -  (1  -  u)£p(N  +  p) 

(  N  -1 

+  (1  -  u)P j-a(p  +  r)  l  (N  -  n)(p  +  n  -  l)c  jU11-  (1  -  u) 

(  n»l 

N 

+  [A  (p  +  r)  -  (h  +  p)]  l  (N  -  n)c  u 

n-0 
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Consider  only  the  terms  within  the  curly  brackets  which  may  be  rewritten 
using  the  fact  (p  +  n  -  l)cnl  =  nc^  as, 


{  }  =  -a (p  +  r)  l  (N  -  n)(p  +  n  -  1)cr  jU 
n=l 


+  £  [N  -  n]un[a(p  +  r)n  +  X (p  +  r)  -  (h  +  p))cn  • 

n=0 


Collecting  the  coefficients  of  u  ,u  ,  . u  separately  and  simplifying, 


N-l  N-l  n 

{  }  =  a(p  +  r)  l  nc  un  +  X(p  +  r)  [  c  u 
n=0  n  n=0 

.  N-l 

-  (h  +  p)  l  (N  -  n)c  u  . 
n=0 


Hence, 


)  I 


(3.12) 


S  (u)  -  (1  -  u)|  p(N  +  p)  +  (1  -  u)" 


»  N-l  _ 

-  (h  +  p)  l  (N  -  n)c  u 
n*=0 


N-l 


(p  +  r)  l  (an  +  A)c  u" 


n=0 


for  all  ^  <  u  <  u^  and  for  all  N  *  1,2,  ...,  N  . 


Fact  1; 

As  u  -  1  ,  SN(u)  -  0  .  Hence,  SN(T)  +  0  as  T  ->  »  for  all  values 
of  N  =■  1 , 2  ,  . . . ,  N  . 


Fact  2: 

From  (3.7), 


<=> 


FJ(T)  <  0 
F^(u)  <  0 


for  all  T  >  T, 

*  N 

for  all  u  >  u„  . 

*  N 
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Thus,  from  (2.24),  it  follows  that 


FjJj(u)  -  -p  +  (1 


(3.13) 


-  u)P  ( 


N 


(h  +  p)  l  cnun  -  (X  +  aN)  (p  +  r)cNuN 
n-0 


<  0 


for  all  u  >  and  for  all  N  -  0,1,2,  ...,  N  . 


Proposition  3.3: 

For  all  N  -  1,2,  . ..,  N  ,  the  function  S^(T)  (as  defined  in  (3.11)) 
(strictly)  decreases  as  T  increases  in  the  range  [TN_^,TN)  • 


Proof: 

To  show  SN(T)  is  strictly  decreasing  in  T  e  lTN_i»TN)  >  ic  is 
sufficient  to  show  that  S^(T)  <  0  for  all  T  e  [T^,^)  ,  equivalently, 
S^(u)  <  0  for  all  u  e  (u^.u^)  since  S^(T)  »  S^(u)a(l  -  u)  .  Now, 


dS  (u) 

SN(U)  "  dS - (N  +  P)P 


>0+l[a 


N-l  .  N-l  , 

n-1  ...  .  r  n-l 


+  (1  -  u)p  a(p  +  r)  l  nc  nun  +  X(p  +  r)  £  nc  u 
1  n-l  "  n-l 


.  N-l 

-  (h  +  p)  l  (N  -  n)nc 
n-l 


n-l] 

»"  J 


I" 


hN-1  N-l 

+  r)  l  nc  un  +  X(p  +  r)  l  c  u 
n-l  n  n-0 


*  N-l  1 

-  (h  +  p)  l  (N  -  n)c  un 

n-0  J 


Taking  (1  -  u)P  as  a  common  factor  and  using  the  fact 

nc  -  (p  +  n  -  l)c  ,  ,  we  get  after  rearranging  the  terms, 
n  n-l 
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N-l 


(u)  *  -(N  +  p)p  +  (1  -  u)p<  a(p  +  r)  £  (p  +  n  -  ljc^nu 


n-l 


n-l 


-  a(p  +  r) 


N-l  1 

T  nc  nun 

n-l  n  J 


N-l 


+  I A(p  +  r)  l  (p  +  n  -  l)c  ,un  1  -  X (p  +  r) 
1  n-l  n-i 


N-l  I 
I  nc  un 

n-l  n  J 


.  N-i  A  N-l 

-  (h  +  p)  l  (N  -  n)  (p  +  n  -  l)c  .u  -  (h  +  p)  [  (N  -  n)nc 
n-l  n_I  n-l 


N-l  N-l 

-  I (p  +  l)a(p  +  r)  l  nc^u  +  (p  +  l)A(p  +  r)  l  c^u" 


->1 
"  J 


n-l 


N-l 


n-0 


+  (p  +  l)(h  +  p)  l  (N  -  n)c  un 
n-0  n 


concentrating  only  on  the  terms  within  the  curly  brackets,  collect  the 
coefficients  of  u^u1,  uN  ^  and  after  simplifying, + 


S^(u)  -  (N  +  p) ^ (1  -  u)' 


(3.14) 


T  *  N-i 

(h  +  p)  l  c  u 

L  n-0  n 

•  1- 

N-l 


-  [A  +  a(N  -  1) ]  (p  +  r)cN_xj  -  p 
for  all  u^^  m  u  <  ^  ant*  ^or  ^  “  1,2,  ...,  N 


Since  (3.13)  is  true  for  all  N  -  0,1,2,  ...,  N  ,  it  is  true  in  particular 
for  N  -  N  -  1  e  [0,N]  also.  Thus,  from  (3.13)  at  N  -  (N  -  1)  , 


^This  process  is  little  bit  involved.  We  suggest  writing  the  coefficients 

un  separately  for  all  n  -  0,1,2,  ...,  N-l  and  then  cancel  the  terms 
using  the  fact  p  -  A/a  . 
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A  -  N-l 

F^_l<u)  -  -p  +  (1  -  u)P  (h  +  p)  £  cru 

L  n=0 

-  [X  +  a(N  -  1) ]  (p  +  r)c^,  ^  <  0  for  all  u 

Thus, 

(3.15)  S^(u)  <  0  for  all  u  >  u^_1  . 


In  particular, 

S^(u)  <0  for  u  e  [uN_i,uN]  . 

Since  a  >  0  and  (1  -  u)  >0  for  all  T  >  T  (>  0)  , 

*  o 

SJ<T)  <  0  for  all  T  e  [T^.y 

and  for  all  N  *  1,2,  ...,  N  . 

Proposition  3.4: 

For  all  N  -  1,2 . (N-l)  ,  the  function  SN(T)  has  a  positive  jump 

«  T  ■  tn  <  ’  • 


Proof ; 


It  is  sufficient  to  show  S^(u)  has  a  positive  jump  at  u  *  uN  <  1 


4W  ■  WV  -  W  • 


By  (3.12), 


■  V 
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-  p(l  “  V  +  (1  -  Ujj) 

+  X(p  +  r) 


P+1 


a(p  +  r) 


N  N-l 

r  nr  n 

)  nc  u„  -  >  nc  u., 

n  N  n  w 


n-0 


n-0 


N  N-l 

r  n  r  n 

.  I  CnUN  -  l  CnUN 

Ln-0  n=0 


N  N-l 

-  <h  +  P)|  l  (N  +  1  -  n)cnuJJ  -  l  (N  -  n)cnuJJ 
Ln-0  n-0 


Cancelling  some  of  the  terns  after  expansion. 


AS, 


N(V  "  (1  ”  'V 


-  (h  + 


\ 

-  N  nl 

p)  y  c  uJ 
n-0  nNJ 


(X  +  aN)  (p  +  r)cNuN 


By  (3.13), 


(1  -  u)' 


r  -  n 

(h  +  p)  l  c  un  - 

L  n-0  n 


(X  +  aN) (p  +  r)cNu 


N 


-  p  <  0 


for  u  >  and  for  all  N  -  0,1,2,  ...»  N  . 


Thus, 


4VV 


>  0 


for  all  N  -  1,2,  . . . 
and  u^  <  1  . 


Proposition  3.5: 


For  N  -  1,2,  ...,  N  (finite), 


(i)  tt^(T)  can  have  at  most  one  change  in  sign  in  the  range 

T..  .  <  T  <  Tx,  . 

N-l  •  N 

a 

(ii)  In  particular,  if  it^(T)  <  0  for  some  T  z  [1^,1  )  * 
•rr^ (T)  is  strictly  concave  thereafter  up  to  TN  . 


,  (N-l) 


then 


Proof : 


(i)  From  (3.10), 


■;<t> 


“Xp  -  pp  +  Sn(T) 


By  Proposition  3.3,  S^(T)  strictly  decreases  in  the  range 

T.,  ,  <  T  <  T.t  .  Hence,  once  the  numerator  becomes  negative  in 
N-l  *  N 

that  interval,  it  stays  negative  in  that  interval,  since  p  >  0 


and  (-Xp  -  pp)  is  a  constant. 

A  A 

(ii)  Assume  tt  '  (T)  <  0  for  T  e  [T  ,  ,T  ) 
N  m  N-l  N 


*”(T)  -  ,;•(«)(!  -  u)a  , 


since  tt'(T)  =  tt'(u)  .  From  (3.10), 
N  N 


"n(u) 


-Xp  -  pp  +  SN(u) 
1  -  u 


Hence, 


*i’ (T> 


u)S^(u)  +  (-Xp  -  pp  +  SN(u))](l 
(1  -  u)2 


for  all  T  .  <  T  <  T  . 

N-l  “  N 


By  (3.15),  S^(u)  <  0  for  u  2  •  By  assumption 


"n(t)  ^  0 


-Xp  -  pp  +  SM(u)  <  0  . 


Hence, 
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ir’’(T)  <0  for  all  T  <  T  <  T„ 

N  ■  N 

A 

*>  7n(t)  1s  strictly  concave  in  the  range  T  <  T  <  . 

Corollary  1: 

For  N  •  1,2,  ...,  N  ,  tt^T)  attains  its  maximum  value  for  a 

T  e  ^N-l’V  ' 

Proof : 

Follows  directly  from  Proposition  (3.5). 

Proposition  3.6; 

The  expression  (T)  (refer  to  (3.11))  is  strictly  positive 
N  -  1,2 . N  and  T  e  [T^.y  . 

Proof; 

From  (3.11)  and  (3.12), 

VT)  ‘  Vu) 

r.  r  n-i 

-  (1  -  u)  p(N  +  p)  +  (1  -  u)P  (p  +  r)  J  (an  +  A)c  u11 

L  L  n-0  n 

-  (h  +  p)  l  (N  -  n)c  un|l 

n-0  JJ 

for  all  uN_^  <  u  <  uN  and  N  -  1,2,  ...,  N 


By  expanding  the  terms,  the  following  identity  may  be  verified 


N-l  _  N-l  N-2 

I  (N  -  n)cnu  =  l  cnun  +  l  c/  +  ... 
n-0  n-0  n-0 


2  1 

+  y  c  un  +  y  c  un  +  c  u°  . 

n-0  n  n-0  n  0 


unique 


for  all 


Thus,  S^(u)  may  bo  rewritten  using  the  above  Identity  as, 


SN(u)  -  (1 


(  N-l 

-  u)  p(N  4  e)  +  (1  -  u)DUp  ■+  r)  l  (an  +  X)c  u° 

(  n-0 


-  N-l  /  m  \) ' 

-  (h  +  p)  l  (  I  cA 

m-0  \n=0  /)J 


for  u„  -  <  u  <  u,,  and  N  «  1,2,  ...,  N  . 
N-l  ■  N 


After  rearranging  the  terms, 


p  4-  (1  -  u)°  (p  +  r)(am  +  A)c  um 

tn 


(3.16) 


Sj^Cu)  -  (1  -  u)£pp  +  l  p  +  (1  -  u)°  |(p  +  r)  (am  • 

'  <h  +  P)  Jo 

for  all  uJJ_1  <  u  <  u^  and  N  =  1,2,  ...,  N  . 


Since  u  -  1  -  e 
n 


for  ail  n  »  0,1,2,  ...,  N-l  ,  and 


T_  <  T,  <  ...  <  T 


'  (refer  to  Corollary  3),  u_  <  u.  <  . . 
N-l 


using  (3.13),  for  all  uN_^  <  u  <  uN  , 


Hence, 


(3.17)  V  p  4-  (1  -  u)P  (p  4-  r)(am  4-  A)c  um  -  (h  +  p) 


.  m  1 

P)  I  c  u 
n-0  ” 


>  0  . 


Using  (3.17)  and  the  fact  op  >  0  and  (1  -  u)  >  0  , 


SN(u)  >  0 


for  all  u.,  <  u  <  u 

N-l  -  N 


and  N  -  1,2,  . . . ,  N 


Sn(T)  >  0 


(3.18) 


for  all  Tn_x  <  T  <  Tn 


and  N  -  1,2,  ...,  N  . 
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Proposition  3.7: 


There  exists  a  finite  lower  bound  >  To  for  the  review  period  such 


k  +  C(y  <T),T)  -  G(0,T) 


for  T  “  T 
>  £ 


for  finite  set-up  cost  k  . 


Proof: 


By  (2.13),  yQ(T)  -  0  for  0  <  T  <  Tq 


k  +  G(yo(T),T)  >  G(0,T) 


for  0  <  T  <  T  , 

*  o  ’ 


since  G(y  (T)  ,T)  =  G(0,T)  .  Hence,  T  t  [0,T  )  .  So,  consider  only 

O  jL  0 

positive  order  levels,  i.e.,  y  (T)  >  1  (=>T  >  T  )  ,  At  T  *  T  , 

o  *  *  o  o 

G(yo(T),T)  -  G(l,To)  -  G(0,To)  since  £,G(0,To)  is  zero.  (Refer  to 


Section  2.13,  Corollary  3.) 

dG(y  (T),T) 
Let  G'(yo(T),T)  - - - 


.  From  (3.1), 


(3.19) 


G'(yo(T),T) 


-  *’(T) 


for  all  T  <  T  <  ®  . 

o  * 


Note: 


In  (3.19),  only  right-hand  derivatives  exist  at  Tn,T. ,  ...,  T 

N-l 

since  y  (T)  lias  finite  saltus  at  these  points  and  is  right  continuous. 
Substituting  the  value  of  ’r'(T)  from  (3.10)  and  using  the  fact 


y„<T>  -  N  for  T„_1  i  I  <  T„  . 


A  (p  +  r  -*•  p/~\)  -  S  (T) 

G’(y  (T),T)  -  C'(N,T)  - - - - 


for  j  <  T  <‘T,3  and  N  •  1,2,  ....  N 


0. 20) 


Hence, 


Xl  'XI  *XI 
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(3.22)  G'(0,T)  -  G'(yo(T),T)  >  0  for  all  T  e  [Tq,») 

since  yQ(T)  ■  N  for  <  T  <  and  N  ■  1,2,  ...,  N  .  Since 

G(yo(T),T)  and  G(0,T)  have  the  same  value  at  T  ■  Tq  and  the  rate  of 
increase  of  G(0,T)  is  strictly  greater  than  that  of  G(yQ(T),T)  (by 

(3.22) ),  the  difference  [G(0,T)  -  G(yo(T),T)]  strictly  increases  in  T  . 

Once  the  difference  equals  the  finite  set-up  cost  k  at  T  ■  T£  ,  it  will 

continue  to  increase  (strictly)  for  all  T  >  T.  (>  T  )  ,  thus  proving  the 

X  O 

proposition. 

Comment  1: 

A  positive  inventory  level  is  required  to  remain  in  business,  i.e., 
yQ(T)  >  1  .  Thus,  the  value  of  the  lower  bound  on  the  review  period  T£  is 
necessary  since  y  (T)  £  0  for  T  e  [0,T£]  •  Let  denote  yQ(T£)  • 

Hence,  an  optimal  review  period  that  exceeds  is  sought  so  as  to  get  an 

order  level  at  least  equal  to  N£(>  1)  . 

Note; 

(3.23)  y0(T)  ■>  N£  for  T£  <  T  <  . 

3.6  Algorithm  to  Compute  Optimal  Review  Period 

It  was  shown  in  Section  3.2  that  n(T)  tends  to  -*  as  T  ■*  *  .  So, 

n(T)  reaches  its  maximum  value  at  some  T  <  ®  .  Corollary  1  to 

Proposition  3.5  implies,  for  N  -  1,2,  ...,  N  ,  there  exists  a  unique  T 

for  every  interval  [T.,  ,,T„)  where  %,(T)  attains  its  maximum  in  this 

N-l  N  N 

interval.  Since  the  number  of  intervals  are  finite,  the  problem  reduces  to 
just  finding  the  maximum  value  of  tt^(T)  for  each  N  *  ••••  ®  and 

choosing  the  best  of  all  the  maximuos.  (Refer  to  Comment  1  and  (3.23)  at 
the  end  of  Section  3.5.) 
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Let  T  denote  the  point  in  the  int  jrval  [T,  .  ,T..)  where  it  (T) 

N  N-l  N  N 

attains  its  maximum  value.  Any  one  of  the  following  three  situations  may 
arise : 


‘l>  ;  0 


Is  T 


By  Proposition  3.5,  tt^(T)  decreases  in  this  interval.  Hence, 
rN-l  * 


(ii)  0  and  "n(tn)  ■*  0  • 


Again  by  Proposition  (3.5),  ^(T)  ■>  0  for  all  T  c  [T^,  ^,T^,)  •  Hence, 

tt.,(T)  strictly  increases  in  this  interval  and  T„  *  Tv.  . 

N  N  N 


(ill)  0  Jnt*  *J<V  <  0  • 


By  Proposition  (3.5)  and  Corollary  1,  there  exists  a  unique  T,  where 

tt'(T%.)  *  0  •  Hence,  the  maximum  value  is  attained  at  Tv  . 

N  N  ft 

A  step  by  step  algorithm  to  compute  the  optimal  review  period  is  given 
which  can  be  programmed  directly  on  a  digital  computer. 


Algorithm: 


Denote  the  optimal  review  period  by  T  and  c he  optimum  value  of  the 
order  level  optimized  net  return  function  by  tt  .  Let  N  ■  N  . 


Initializing  ften:  (Refer  to  Comment  1  and  (3.23)  at  the  end  of  Section  3.5.) 


Compute  tt *  (T  ;  .  If  ”’(1^)  s  0  ,  then  set 


m(T  )  ,  T  *  T„  and 

L  K 


go  to  Step  2. 


If  t'(T\)  >  0  ,  there  are  two  possibilities: 


(i)  m’(T;)  >  o  ,  j  >  0  .  T1 


Then  set 


W)  • 1  ’  \ 
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go  to  Step  2. 

(ii)  ir’(Tt)  >  0  ,  n^TN  j  <  0  .  Then,  seek  a  T  e  ^T£,Tn  J  such 

that  it*  (T)  is  zero.  Set  n  -  tt(T)  ,  T  -  T  and  go  to  Step  2. 

Step  1; 

Compute  .  If  <  0  ,  go  to  Step  2.  Otherwise,  go 

to  Step  3. 

Step  2: 

Increase  K  by  one.  If  N  <  N  ,  go  to  Step  1.  Otherwise,  go  to 
Step  3. 

Step  3: 

Compute  n '  (T^)  .  If  ir^CT^)  <  0  ,  go  to  Step  4.  Otherwise,  set 
it  ■  ^TN>  *  ^  *  TN  and  8°  C°  StCp 

Step  4: 

Compute  T  in  the  range  T^_^  <  T  <  IN  such  that  1T^(T)  *  0  • 

(Since  in  this  case  n'( T  )  >  0  and  <  0  ,  the  bisection 

N  w-1  N  N 

procedure  described  in  Chapter  2  may  be  used  to  find  T  .)  Compute  r  (T) 
*  ^  » 

If  n^(T)  <  v  t  go  to  Step  2.  Otherwise,  set  n  ■  t^(T)  ,  T  ■  T  and  go 

to  Step  2. 

Step  5: 

Terminate.  The  current  value  of  T  is  the  optimal  review  period  and 
ir  is  the  maximum  value  of  the  profit  function. 

Finiteness  of  Algorithm 

By  Proposition  2.6,  the  order  level  is  bounded  as  i  -  «  (i.e.,  N 

finite).  Hence,  only  a  finite  number  of  intervals  [T^^.T^)  for  ^ 


is 


8? 


ranging  from  •••»  N  aro  to  be  examined  for  the  maximum  value  of 

ti(T)  .  Hence,  the  algoritlim  terminates  in  a  finite  number  of  stops. 

3.7  Case  When  "Lost  Sales"  Costs  are  Ignored 

In  Section  3.0,  it  was  mentioned  that  the  cost  of  "lost  sales"  may  be 
ignored  by  rnakiug  p  and  p  zero.  This  may  simplify  the  cost  and  revenue 
expressions  considerably.  Since  p  and  p  are  assumed  to  be  nonnegative 
throughout  the  analysis  so  far,  all  the  results  of  optimal  order  policy  and 
optimal  review  period  hold  for  this  case  also.  An  example  given  in 

Figure  3.1  illustrates  the  graph  of  F^(T)  against  T  for  N  *  0,1,2 . 6 

when  p  and  p  are  zero  and  h  is  strictly  positive. 

3.8  Treatment  of  Back  Order  Regime 

It  was  pointed  out  in  Section  3.0  that  one  can  follow  the  policy  of 
filling  up  the  back  orders  at  the  end  of  period.  Since  it  is  a  single-period 
problem,  the  revenue  from  filling  the  back  orders  should  be  included  in  the 
revenue  expression.  This  means  an  additional  revenue  equals  (r  -  c)  times 
the  amount  back  ordered,  neglecting  the  set-up  cost  k  for  reorder  ing. 

Thus,  the  net  return  function  will  be 

v(y,T)  “  nn(T)  -  k6(y  -  0)  -  G(y,T) 

m 

+  (r  -  c)  l  (n  -  y)P  (T)  . 
n**y+l 

The  above  expression  may  be  directly  obtained  from  (2.4)  by  simply  replacing 
the  parameter  "p"  by  "(p  +  c  -  r)."  Thus,  all  the  results  of  optimal 
order  policy  and  its  properties,  given  in  Chapter  2,  follow  identically  with 
the  replacement  of  "p"  by  "(p  +  c  -  r)." 

Optimal  Review  Period 

Though  the  analysis  of  Chapter  2  went  without  any  hitch,  it  does  not 


FICURE  3.1 


happen  so,  for  finding  an  optimal  review  period.  The  place  where  it  differs 


is  in  the  limit  of  the  order  level  optimized  revenue  function  as  T  ■*  »  . 
In  other  words, 

lim  7T (T)  *>  -«  is  not  true  always. 

T-x» 


Referring  to  Section  3.3,  and  replacing  "p"  by  "p  +  c  -  r"  ,  the 
expression  for  tt (T)  may  be  written  as, 

»<T)  -  z  (Q:  +  Q2  +  Q3  +  Q4) 

P 


where 


(3.11),  with  the  substitution  p  ■  (p  +  c  -  r)  .  Thus, 


(3.24) 


where 


•i(T) 


X[r  -  (p  +  c)]  -  pp  +  Sn(T) 


Sn(T) 


(3.25) 


P  (N  +  p)p  -  (h  +  p)  l  (N  -  n)P  (T) 
1  n=0  n 


(p  +  c)  l  [N  -  n][(p  +  n  -  l)paP  , (T) 
n“0 


for  all  T„  .  <  T  <  T.,  . 

N-l  =  N 


Comparing  (3.11)  and  (3.24),  it  is  clear  that  Propositions  3.3  thru  3.7 
hold  for  this  case  too.  Also,  (3.18)  holds,  i.e., 


sN(T)  >  o 


(3.26) 


for  all  Tm  .  <  T  <  TM 
N-l  *  N 


and  N  *  1,2,  . . . ,  N  . 


If  X[r  -  (p  +  c)]  -  pp  <  0  ,  then  the  results  of  Section  3.5  hold  and 
the  algorithm  to  compute  the  optimal  review  period,  as  given  in  Section  3.6, 
may  be  applied  without  modification.  On  the  other  hand,  if 
X[r  -  (p  +  c)]  -  pp  >  0  ,  then 


irJ(T)  >  0 


for  all  T  e  [T  ,»)  . 


(See  (3.24)  and  (3.26).)  Hence,  n (T)  increases  indefinitely  as  T  increases 
and  the  optimal  review  period  may  be  made  as  large  as  possible.  This  case 
may  not  happen  in  a  number  of  systems  as  "c"  may  be  as  high  as  "r,"  if 


the  product  is  bought  from  a  competitor  so  as  to  satisfy  the  customers,  and 
the  back  order  cost  "p"  in  general  is  positive. 
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Since  seasonal  goods  have  a  finite  upper  bound  for  the  length  of  the 
review  period  (i.e.,  the.  length  of  the  season),  even  if 
X(r  -  (p  +  c)]  -  pp  >  0  ,  the  optimal  review  period  can  only  be  made  as 
large  as  the  length  of  the  season. 
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CHAPTER  4 

ESTIMATION  OF  PARAMETERS 

4.1  The  Estimation  Problem 

As  discussed  in  Chapter  1,  the  probability  distribution  of  demand  is 
given  by 


Pn(T) 


r(P  +  n) 
r(p)r(n  +  i) 


,  -aT.p,,  -aT.n 
(e  )  (1  -  e  ) 


where  Pn(T)  *  probability  of  n  demands  in  [0,T]  and  p  ■  A/a  >  0  . 

Once  the  parameters  A  (the  constant  demand  rate)  and  a  (the  unit 
contagion  rate)  are  determined,  Pr(T)  is  known  for  given  T  ,  and  the 
results  of  Chapter  2  and  3  may  be  applied  to  find  optimal  order  policy  and 
review  period.  Thus,  we  are  left  with  estimating  the  values  of  a  and 
A  using  the  knowledge  of  the  demands  in  the  previous  periods.  Of  course, 
this  poses  a  problem  for  the  first  period  as  there  is  no  prior  knowledge. 
This  can  be  overcome  if  a  sample  survey  is  done  where  the  product  was 
introduced  in  a  smaller  scale  in  a  sample  area.  As  a  matter  of  fact,  the 
practice  of  conducting  a  pilot  study  is  prevelant  in  many  cases.  For 
example,  a  cereal  manufacturer  usually  introduces  the  product  in  a  small 
sample  area  before  it  is  introduced  nationwide.  Using  the  same  estimation 
procedures  described  in  later  sections,  a  good  estimate  of  the  initial 
values  of  a  and  A  may  be  obtained. 

A  number  of  estimation  procedures  are  available  in  literature  like 
the  method  of  moments,  method  of  frequencies,  maximum  likelihood  procedure, 
modified  chi-square  method  and  so  on.  We  shall  essentially  describe  two 
methods  and  exhibit  their  estimators.  The  first  method  will  be  the  direct 
and  simple  method  of  moments  which  may  not  be  efficient.  But  the  second 


1 
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method  will  be  the  maximum  likelihood  estimate  which  we  will  show  for  large 
sample  size,  that  it  is  asymptotical  ly  unbiassed,  minimum  variance  and 
efficient  estimator  of  the  parameters. 

A,  2  Method  of  Moments 

This  is  the  oldest  general  method  proposed  for  estimating  the  values 
of  the  parameters  of  a  distribution  by  means  of  a  set  of  sample  values. 

This  method  consists  in  equating  a  convenient  number  of  the  sample  moments 
to  the  corresponding  moments  of  the  distribution,  which  are  functions  of 
the  unknown  parameters.  By  considering  as  many  moments  as  there  are 
parameters  to  be  estimated  and  solving  the  resulting  equation  with  respect 
to  the  parameters ,  estimates  of  the  latter  are  obtained.  This  method  often 
leads  to  comparatively  simple  calculations  in  practice. 

Assume  we  have  N  sample  areas  and  random  variables  X^.X^, 
denote  the  number  of  demands  in  [0,T]  from  the  ith  sample.  X^'s 
are  independent,  identically  distributed  and  have  a  contagious  dist ribut 
such  that 

«.n  pv>(t)  ■  rdmtfn  (C'“T)°  (1 *  e"T)X  £or  311  x '  °'1'2- 

Define: 

N 

l  h 

i*l 

Sample  Mean  *  X  =  - 

2 

Sample  Variance  =*  s  * 


l  (x,  -  X) 


i-1 


j 


We  know  (from  Chapter  1)  that  the  theoretical  moments  are  given  by, 


/•.  -aT. 
p (1  -  e  ) 
m(T)  -  expectation  -  — 3^, - 


,,  -aT. 

V(T)  «*  variance  =  ^ 


Equating  the  theoretical  and  sample  moments, 


N 

l  Xi 


(4.2) 


o(T) 


i«l 


(4.3) 


N  - 

l  (X  -  x)2 

_  i-1  2 

V(T)  .  - - - ■  s 


Dividing  (4.2)  by  (4.3)  we  get,  a  the  estimator  of  a 


-aT  X 

e  -  ,  i.e. , 

s 


°  "  T  Kl 


N  ,, 

l  (X  -  X) 
i-1 

N 

l  h 

i-1 


Substituting  a  in  (A. 2),  we  get  the  estimator  for  X  as 


(4.5) 


X  = 


aX 


sfc  -  X 


as , 


Thus,  a  and  X  are  the  estimators  for  a  and  X  and  the  actual 
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estimates  can  be  found  by  replacing  the  random  variables  . . . ,  X^ 

by  observed  values  ,  x^  in  expressions  (4.  4)  and  (4.5) 

respectively. 

4.3  The  Method  of  Maximum  Likelihood^ 

From  a  theoretical  point  of  view  this  is  the  most  important  general 
method  of  estimation. 

Using  the  same  notations  of  Section  4.2,  we  define  the  likelihood 
function  L  of  the  sample  of  n  values  from  a  population  of  the  discrete 
type  by  the  relation 

(4.6)  L(X. ,  ....  X  ;a,X)  -  P  (a,X)  ...  p  (a,X)  . 

1  n  x,  x 

1  n 


When  the  sample  values  of  X.,  ....  X  are  given,  the  likelihood  function 

i  n 

L  becomes  a  function  of  the  two  variables  a  and  X  .  The  method  of 
maximum  likelihood  now  consists  in  choosing,  as  an  estimate  of  the  unknown 
pair  of  values  (a,X)  ,  the  particular  value  (a,X)  that  maximizes  the 
likelihood  function.  Since  log  L  attains  its  maximum  for  the  same  value 
of  (ct,X)  as  L  ,  we  thus  have  to  solve  the  likelihood  equation 


(4.7) 


61ogL 

6a 


0 


(4.8) 


6X 


with  respect  to  a  and  X  .  We  shall  disregard  any  root  of  the. form 
a  -  constant  and/or  X  »  constant,  thus  counting  as  a  solution  only  a 
root  which  effectively  depends  on  the  sample  values  . X^  .  Any 

*f»  ^ 

For  a  detailed  discussion  the  reader  is  referred  to  Cramer  [2]. 
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solution  of  the  likelihood  equation  will  then  be  called  a  maximum  likelihood 
estimate  of  (a, A)  . 

The  importance  of  the  method  is  clear  by  the  following  two  propositions: 

*  * 

(i)  If  a  pair  of  joint  efficient  estimates  a  and  A  exists,  the 

* 

likelihood  equations  will  have  the  unique  solution  a  -  a 

A 

and  A  -  A 

(ii)  Any  solution  of  (4.7)  and  (4.8)  is  a  sufficient  estimate  of 
(a, A)  if  a  sufficient  estimate  exists.  From  (4.6), 

n 

log  L  -  l  log  P  (T)  , 

6  i-1  6  Xi 

and  from  (4.1) 

log  Px  (T)  -  log  T (p  +  xi)  -  AT  +  x^  log  (1  -  e  aT) 

-  log  T(p)  -  log  r(xt  +  1)  . 


Thus, 


-  —  dT 

log  L  =  £  log  r(p  +  x  )  -  nAT  +  nx  log  (1  -  e  ) 

e  i-1 


n 

-  n  log  r(p)  -  l  log  r(x  +  i)  . 
i-1 


Define 


4<(y) 


A  d  log  r(v)  A  r'(v) 
(Ay  "  r  (y) 


(4.10) 
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ij*(y)  is  known  in  literature  as  digamma  function  and  tables  are  available 
in  Fletcher,  et  al  (refer  Appendix  III  on  Gamma  functions).  Setting 

"  0  and  using  (4.10),  we  get 


(4.11)  “  l  *(P  +  x  )  -  nT  -  J  *(p)  -  0  . 

i-1 


Similarly, 


61ogL  1 
6a  2 


<KP  +  xt)  + 


nXTe 


-aT 


0  . 


Using  (4.11),  this  reduces  to 


A 

a 


» 


{ 


i*6t  i 


(4.12)  X(eaT  -  1)  ■  ax  . 

Now  we  are  left  with  two  simultaneous  equations  in  two  unknowns.  Of  course, 
a  direct  approach  is  to  solve  (4.11)  and  (4.12)  iteratively  for  a  pair 
(a,X)  .  This  should  not  present  any  problem  once  it  is  committed  to  a 
digital  computer. 

4.4  Efficiency  of  the  Maximum  Likelihood  Estimator 

We  shall  show  a  few  properties  about  the  contagious  distribution  given 
in  (4.1)  which  are  necessary  to  show  that  the  pair  of  values  (a, A)  obtained 
from  (4.11)  and  (4.12)  is  asymptotically  unbiassed,  efficient  estimate  of 
(a, A)  as  the  number  of  samples  is  large. 
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PI: 

The  first,  second  and  third  derivatives  of  log  Pn(T)  with  respect 
to  a  and  \  exist. 

Proof : 

The  truth  of  this  statement  follows  from  the  fact  that  Gamma  function 
has  continuous  derivatives  of  all  orders  (refer  Appendix  III), 

P2: 

The  contagious  distribution  { (T) }  has  finite  moments  of  all  orders, 
for  given  T  . 

Proof : 

For  given  T  ,  the  generating  function  of  pn(T)  is  given  by, 


Similarly, 


9‘, 


G"(z>  - 

<1  -  qZ) 


.2 


and  G"(l)  «  £  n(n  -  1)P  (T)  *  ^.lP— t.-U.q..  #  Thus,  in  general, 


n-0 


Gn(l)  -  Lt  — ^  -  l  n(n  -  1)  ...  2. 1*P  (T) 
Z-l  dZ  n-0 

_n 

m  P(P  +  1)  ...  (o  +  n  -  l)q 

-n 

P 


Thus,  P  (T)  has  finite  moments  of  all  orders  for  Te  [0,®)  . 


P3: 


l  r(p  +  n)P  (T)  is  finite  and  positive. 
n-0  n 


Proof : 

It  is  clear  that  the  summation  is  positive  as  T(p  +  n)  >  0  for  all 
n  -  0,1,2,  ...  and  p  >  0  .  Finiteness  follows  from  P2. 


The  reader  can  verify  himself  that  Cramer's  conditions  for  asymptotic 
efficiency  of  maximum  likelihood  estimators  follow  directly  from 
Properties  PI,  P2  and  P3.  Thus,  the  maximum  likelihood  estimator  (A, a) 
given  by  (4.11)  and  (4.12)  is  asymptotically  unbiassed  efficient  estimator. 
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APPENDIX  I 

SOLUTION  OK  ORDINARY  DIFFERENTIAL  EQUATION 


Reference: 

Granville,  W.A. ,  P.F.  Smith  a  ml  W.R.  Long ley,  ELEMENTS  OF  THE 
DIFFERENTIAL  AND  INTEGRAL  CALCULUS,  Ginn  and  Company,  Snn  Francisco,  n.  380, 
(1934). 


Here  we  give  n  method  to  solve  ordinary  differential  equations 
of  the  form 


(1) 


dx 


+  Py 


Q 


where  P  and  Q  are  functions  of  x  alone  or  constants. 
To  integrate  (1),  let 


(2) 


y  *  uz 


where  z  and  u  are  functions  of  x  to  be  determined.  Differentiating  (2), 


(3) 


±L 

dx 


dz 

u  ~  +  z 
dx 


du 

dx 


Substituting  from  (3)  and  (2)  in  (1),  we  get 


u 


dz 

T*  +  z 
dx 


du 

dx 


+  Puz 


Q  , 


i .  e . , 


(4) 


We  now  determine  u  by  integrating 


+  Pu  *  0  , 


in  which  the  variables  x  and  u  are  separable.  Using  the  value  of  u 
thus  obtained,  we  find  z  by  solving 


u  e 

ud^- q  • 


in  which  x  and  z  can  be  separated.  Obviously,  the  values  of  u  and 
z  thus  found  will  satisfy  (A),  and  the  solution  of  (1)  is  then  given  by 


The  differential  Equation  (1.5)  we  had  in  Chapter  1,  is  given  by 


P'(t)  *  -(X  +  an)P  (t)  +  [X  +  a(n  -  1)]P  . (t)  . 

n  n  n-i 


By  putting  y  ■  Pn(t)  ,  x  ■  t  ,  P  *  (X  +  an)  , 


Q  »  X  +  a(n  -  1)  , 


we  will  have  the  same  form  as  in  (1).  Hence,  we  use  the  method  described 


just  now  to  solve  for  P  (t) 

n 


Tables  of  Incomplete  Beta  functions  and  ratios  are  available  in  Pearson 

[11]. 

It  can  be  verified*  using  (1)  and  (2)  when  m  is  an  integer, 

that 


I  (m,n)  *  1  -  (1  -  x)n  l  . 

i-0  '  ' 

The  advantage  of  this  equation  is  that  of  computational  feasibility.  The 
above  series  is  made  up  of  entirely  positive  terms  and  hence  can  be  summed 
quite  accurately,  even  for  fairly  large  values  of  parameters  m  and  n  . 


The  actual  derivation  is  given  in  Harter,  H.L.,  NEW  TABLES  OF  INCOMPLETE 
GAMMA  FUNCTION,  ETC.,  Aerospace  Research  Laboratories,  U.S.  Government 
Printing  Office,  Washington,  D.C.  20402,  (1964). 


A. 5 


APPENDIX  ITT 
THE  GAMMA  FUNCTION 


Reference: 

Cramer,  H.  [2]. 

The  Gamma  function  T(p)  is  defined  for  all  real  p  >  0  by  the 

integral 

00 

(1)  T(p)  -  J  xP  1e  *dx  . 

0 

The  function  is  continuous  and  has  continuous  derivatives  of  all  orders: 

00 

(2)  r(r)(p)-  y  xP_1(log  x)re“Xdx 

0 

for  any  p  >  0  .  When  p  tends  to  0  or  to  4®  ,  T(p)  tends  to  4»  , 
Since  the  second  derivative  is  always  positive,  T(p)  has  one  single 
minimum  in  (0,®)  .  Approximate  calculation  shows  that  the  minimum  is 
situated  in  the  point  pQ  ■  1.4616  ,  where  the  function  assumes  the  value 
r(p  )  -  0.8856  . 

By  a  partial  integration,  we  obtain  from  (1)  for  any  p  >  0 
T(p  +  1)  =  pT(p)  . 

When  p  is  equal  to  a  positive  integer  n  ,  a  repeated  use  of  the  last 
equality  gives,  since  T(l)  *  1  , 


T(n  +  1)  ■  n  ! 


From  (1),  we  further  obtain  the  relation 


f  x»-V“dx  -  ^ 
0 


where  a  >  0  ,  X  >  0  .  From  (2),  when  r  *=  1  , 


r’(P) 


/ 


XP  1(log 


x)e  Xdx 


Define  ii>(p)  *  ^ (p)  is  called  digamma  function.  Table 

dp  i  ip; 

of  digamma  functions  are  available  in  Fletcher,  A.  et  al,  AN  INDEX  OF 
MATHEMATICAL  TABLES  I,  Sect.  14. A,  Addison-Wes ley  Publishing  Co.,  Inc., 
Reading,  Massachusetts,  (1962). 
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